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Abstract 



Two approaches to quantum-classical correspondence are distinguished according to the 
classical dynamical theory with which quantum theory is compared. The first of these, 
Ehrenfest correspondence, defines a dynamical regime in which the quantum expectation 
vàlues foUow approximately a classical trajectory. The second of these, Liouville correspon- 
dence, applies when the quantum probability distributions remain well approximated by a 
density in the classical phase space. The former applies only for narrow states, whereas the 
latter may remain vàlid even for quantum states that have spread to the system size. 

A spin model is adopted for this correspondence study because the quantum state is 
discrete and finite-dimensional, and thus no articifical truncation of the Hilbert space is 
required. The quantum time-evolution is given by a discrete unitary mapping. The cor- 
responding classical model is volume-preserving (non-dissipative) and the time-evolution is 
given by a symplectic map. 

In classically chaotic regimes, the widths of initially narrow quantum states grow, on 
average, exponentially with time, until saturation at the system size. This initial spreading 
rate is well approximated by the classical Lyapunov exponent when the accessible classical 
phase space is predominantly chaotic. Because of the exponential growth rate of the quan- 
tum variance, the Ehrenfest regime is delimited by a break-time that grows logarithmically 
with increasing quantum numbers. 

The small differences between quantum expectation vàlues and corresponding Liouville 
averages of dynamical variables also grow exponentially, initially, if the classical behaviour 
is chaotic. This exponential rate is independent of h and consistently larger than the clas- 
sical Lyapunov exponent by at least a factor of two. Interestingly, this exponential growth 
rate does not continue until these differences approach magnitudes of order the system 
size, but crosses over to power-law growth (or simply saturates), when the differences have 
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reached magnitudes no larger than order h. Because of the exponential growth of these 
quantum-classical difFerences, the Liouville regime is also delimited by a break-time that 
grows logarithmicly with increasing quantum numbers. However, due to the early satura- 
tion of the exponential growth of the differences, the Liouville log break-time applies only 
in a restricted domain. 

After spreading to the system size, the quantum states relax to an approximately time- 
independent state. This is characterised by equilibrium distributions that are subject to 
rapidly oscillating fiuctuations about the coarse-grained classical steady-state. If the acces- 
sible classical phase space is predominantly chaotic, then the equilibrium quantum probabil- 
ity distributions are nearly microcanonical. Otherwise the steady-state quantum probability 
distributions accurately reflect the detalls of the KAM surfaces in the classical phase space. 
The equilibrium regime quantum fiuctuations, which are quantum-classical differences, are 
shown to approach zero as a negative power of the quantum numbers. 

Since the differences between quantum expectation vàlues and individual classical tra- 
jectòries grow exponentially in time to a magnitude that scales with the system size, the 
Ehrenfest regime is delimited by an unrestricted log break-time. The log break-time may be 
short compared to experimental time-scales for some chaotic macroscopic bodies, and, as a 
result, an observable breakdown of Ehrenfest correspondence may arise for some macroscopic 
systems. Although the Liouville regime is also delimited by a log break-time, an observ- 
able breakdown of correspondence between quantum mechanics and Liouville mechanics for 
macroscopic bodies is not likely since their differences, expressed relative to the dimensions 
of the observable, approach zero in the limit of large quantum numbers. 
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Chapter 1 



Introduction 



Though it is widely accepted that classical mechanics should emerge from quantum theory 
in an appropriate limit, this "correspondence principle" meets with severe difficulties when 
the classical dynamics exhibit chaos. Research into the physical conditions required for the 
emergence of classical chaos from the underlying quantum description has become a tòpic 
of active investigation, and has led to a number of differing perspectives. 



Ford and coworkers [H^, 43, 44 1 have considered the criterion of algorithmic complexity, 



which has led them to the conclusion that classical chaos can not emerge from quantum 
theory in the macroscopic limit, and, therefore, that quantum-classical correspondence must 
break down for chaotic motion. While Zurek and Paz |96, ^] and others 99 1 

agree that a breakdown of correspondence does arise, they have argued that the required 
correspondence is restored when the effect of a weak-coupling to an environment, a process 



sometimes called decoherence 5C], is taken into account. Alternatively, Ballentine 



Takahashi |91], and Casati and Chirikov |23] have emphasized that correspondence should 
emerge in the macroscopic limit given an appropriate degree of coarse-graining which may 
be considered to represent the finite resolution of macroscopic measurements. The variety 
of perspectives on this issue is due in part to differing views on which observable properties 
must be taken into consideration in order to demonstrate a formal correspondence in the 
macroscopic limit ||7[ |8|, |ïï|, 56 1. 

The detalls of the quantum-classical transition and the quantum manifestations of clas- 
sical chaos are of considerable interest also for the accurate and practical description of 
mesoscopic systems. The onset of widespread chaos in the classical Hamiltonian model 
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of the highly excited states of hydrogen has been shown to accurately predict the unex- 
pected onset of ionization observed in these states under increasing microwave fields |ï^, 64 ] . 
In the classical model, the microwave field itself provides the non-integrable perturbation 
to the otherwise integrable classical dynamics. Under strong magnètic fields, the energy 
level statistics of these Rydberg àtoms exhibit the characterisitic distribution of spacings 
predicted from random matrix theory precisely when the classical model exhibits chaotic 
behaviour |6C]. 

A better understanding of the quantum-classical transition is required also for practical 
computation in the mesoscopic regime. On the one hand, the exact theoretical description of 
even simple microscopic and mesoscopic systems in quantum theory demands computational 
resources (specifically, classical computational resources!) that increase as a power of the 
quantum numbers and exponentially with the number of interacting subsystems. On the 



other hand, the semiclassical techniques developed by Gutzwiller |53] and others become 
computationally impractical due to the proliferation of (unstable) periòdic orbits in the case 
of chaotic motion. Moreover, the accuracy of the semiclassical approximation is not well 
understood in either the time-domain [^] or the energy domain [^], and remains a tòpic of 
active investigation. The results in this thesis suggest that purely classical techniques may 
be applied to accurately describe the time-domain behaviour of chaotic systems even in the 
case of surprisingly small quantum numbers. 

The theoretical Identification of the characteristic properties of quantum systems with 
chaotic classical counterparts has been an interesting and illuminating tòpic of research in 
its own right. This nascent area of research, sometimes called "quantum chaos," is not a 
unified field but encompasses a wide variety of approaches. In analogy with the manner 
in which classical chaos came to be understood, much of the theoretical progress in quan- 
tum chaos has proceeded through the study of idealized model systems, usually taking the 
form of stroboscopic mappings. Some of this research will be reviewed further below. I 
should also mention here that, whereas a fairly complete understanding of the microscopic 
underpinnings of irreversible thermodynamics has recently emerged from the perspective 
of classical dynamical systems theory very little is known about how these char- 

acteristic properties emerge from the more fundamental quantum mechanical treatment of 
these models. The results presented in this thesis will clarify how these statistical proper- 
ties emerge from the underlying quantum dynamics specifically when the classical model 
exhibits chaotic behaviour. 
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Certainly the most intellectually intriguing aspect of this research pertains to the possi- 
bility of demonstrating that experimentally distinct predictions result from two very different 
interpretations of the quantum state. This pedagogical issue is tied to the fohowing prac- 
tical question about quantum-classical correspondence: which classical theory of mechanics 
actually emerges from the quantum description of macroscopic systems? One alternative is 
Hamilton's equations of motion (Newton's laws); another is the classical Liouville equation 
that describes the dynamics of probability densities with non-vanishing support. The for- 
mer describes the exact, deterministic evolution of single macroscopic systems, whereas the 
latter provides only probabilistic predictions for ensembles of similar ly prepared systems. 

The central focus of the work in this thesis is to demonstrate that, for the chaotic mo- 
tion of some macroscopic systems, Standard quantum mechanics may be unable to reproduce 
Newton's laws of motion over experimentally relevant time-scales. Moreover, contrary to 
the conclusion of a recent argument by Zurek , I will provide evidence that this specific 
breakdown of quantum-classical correspondence is not circumvented by taking into account 
the decoherence effects of a stochastic environment. This striking conclusion does not im- 
ply a failure of the correspondence principle, however, since the statistical predictions of 
Liouville mechanics are shown to be compatible with the predictions of quantum theory for 
such systems. The decisive implications of this scenario for the interpretation of the corre- 
spondence principle, as well as the interpretation of quantum theory itself, will be drawn 
out and related to some of the historical debaté on this tòpic. 

In the remainder of this Introduction I will review the relevant literature on quantum 
chaos and quantum-classical correspondence, before undertaking, in the body of this thesis, 
a detailed comparison of the quantum and classical dynamics of a non-integrable model 
system comprised of two coupled spins. The perspective developed from the results of this 
analysis will generalize and extend much of the earlier work on the dynamical properties of 
quantum chaos and quantum-classical correspondence that has been developed through the 
study of different model systems. 



1.1 Classical Trajectòries and Quantum Expectation Vàlues 

A convenient starting point for this subject follows from considering an old theorem due to 



Ehrenfest |^3[ and the oft-repeated correspondence argument that derives from it [ |65| , [8q , 67, 
38|. For a simple Hamiltonian system of the form H = /2m + V{q), the time-dependence 



CHAPTER 1. INTRODUCTION 



4 



of the expectation vàlues {q{t)) = {^p{t)\q\■^p{t)) and {p{t)) = {^p{t)\p\■^p(t)), in an arbitrary 
quantum state ip{t), is prescribed by the differential equations, 

diqjt)) ^ (p(í)) 
dt m 

= (f(<rt). (1.1) 

These equations do not form a closed set since the time-dependence of the force function, 

{F{q)) = {m\F{q)\m) ^ Fm (l-2) 

has not been specified. However, if the operator F{q) = dV/dq can be expanded in a Taylor 
series about the position centroid of the quantum state, {q{t)), then, 



dq^ 



+ ... (1.3) 

g={g) 



where (Ag)^ = (g^) — (g)^ is the state variance. When the quadratic and higher-order terms 
of this expansion are negligible, the system of differential equations ( |1.1| ) becomes closed and 
the expectation vàlues of the dynamical variables reduce to Hamilton's canonical equations 
of motion, 

dqcjt) _ Pcjt) 
dt m 

^ = {F,pJ = F(ç,), (1.4) 

In the above, {•, •} is a Poisson bracket and qc{t) and Pcit) are classical dynamical variables. 

At this point it is natural to ask: under what conditions may we expect this approxi- 
mation to apply? Clearly, if the potential V{q) is either linear or quadratic in g, as for a 
harmònic oscillator, then the quantum equations are identical to the classical ones. More 
generally, this correspondence follows provided the state variance is narrow compared to 
significant non-linear variation in the prescribed force, 

(Ag)^ « - -í^, (1.5) 



and higher-order contributions to the expansion (|l.3| ) are also negligible. The dynamical 
regime in which (|1.5D holds is sometimes called the Ehrenfest regime 0, ^ . 
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A common textbook argument |65, 68] maintains that the approximation ( |1.5| ) should 
be vahd, quite generally, for macroscopic systems. This arguments holds that, for a macro- 
scopic body, the width of a well-locaUsed state may be extremely narrow compared to char- 
acteristic variation in the macroscopic potential (or some other relevant system dimension). 
Consequently, the approximation, 

(Hq)) = Fm, (1.6) 

should hold for a genèric macroscopic system described by such a locahsed state. The 
important conclusion is that the differential equations governing the time-evolution of the 



expectation vàlues {lA) then accurately approximate Hamilton's equations of motion ( |1.4[) . 
Provided it is possible to match the classical initial conditions, quantum theory can then 
describe the deterministic motion of a single macroscopic body by associating expectation 
vàlues with the classical dynamical variables. The requirement on the initial conditions is 
simply that the quantum centroids must be set equal to the initial phase space coordinates 
of the classical trajectory; that is, {q{0)) = ^c(O) and (p(0)) = Pc{0)- 

These observations may be taken as means of justifying the "correspondence principle", 
that is, the view that the predictions of quantum mechanics must coincide with those of 



classical mechanics in cases where the latter is known to be vàlid 37 1. In the present 
context, "classical mechanics" is taken to refer specifically to Hamilton's equations of motion. 
I will call this interpretation of the correspondence principle the Ehrenfest correspondence 
principle. This conjecture may be stated heuristically as, 

K'í-'J^»- 1 for l^oo, (1,7) 

where is a characteristic action and h is Planck's constant. I will provide a more useful 
statement of this conjecture in Chapter 5. 

Of course, the ratio J /h is actually finite for real physical systems. Moreover, initially 
localised quantum states will generally diffuse with time, and the accuracy of the approx- 
imation ( |1.6| ) eventually breaks down. Therefore, to confirm the validity of the conjecture 
( |1.7D as a principle of correspondence, it is important to determine the time-scale on which 
significant deviations between the quantum and classical predictions arise. Let Xc denote 
the generalised classical coordinates, e.g. Xc = (qcPc), and (x) denote the corresponding 
quantum expectation vàlues. Then, the magnitude of their difference 

6(í) = Kx)-Xe| (1.8) 
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is a time-dependent measure of the degree of correspondence. Given some criterion defining 
adequate correspondence, i.e., a threshold X) one may define |25, 57| a time-scale ÍEhr over 
which Ehrenfest correspondence holds from the condition, 



e(í) <x for í < tEhr- 



(l.í 



This threshold-dependent condition delineates the Ehrenfest regime more precisely than the 
approximation (|1.5| ). This condition diíFers shghtly from the "narrow state" condition ( |1.5D , 
though both definitions will coincide, approximately, as a result of Eqs. ( |1.3| ) and if 
the tolerance threshold x ^or the position is set at the characteristic system size, L, given 



by the RHS of (1.5) and the exphcit degree of correspondence for the momentum is ignored. 

The time-scale of Ehrenfest correspondence will generally depend on the detalls of the 
potential and other system parameters, but it is worthwhile to get some sense of the duration 
of this regime by considering the simplest case of a free-particle. For a particle of mass m, 
starting from a Gaussian state with initial Standard deviation Ag(0) = üq and momentum 
hko, the time-dependent probability distribution is given by, 



P{q,t) = \i,{q,t)Y 



■ exp 



2a2(l + í2/^2 



:i.io) 



/2^(1 + íVt2)1/2 

The width of this Gaussian doubles on the time-scale r = 2ma1/h. The time it takes 
this width to grow to a size Aq >> Oq is given by 

Aq Oo m 



t 



n 



(1.11) 



Therefore, in the case of a quantum state describing a macroscopic object, such as a baseball, 
initially localised to within 1 Angstrom, the time it takes this state to diffuse to a width of 
1 cm is of order 10^^ years [5^, 65]. This is much longer than the estimated lifetime of the 



universe! 

Thus it appears plausible to conclude that quantum expectation vàlues can adequately 
describe the deterministic motion of individual macroscopic objects over time-scales much 
longer than the duration of any conceivable experiment. In the following section I will ex- 
plain how the presence of chaos in the classical dynamics challenges this simple picture, as 



well as the adequacy of the conjecture (1.7) as a principle of quantum-classical correspon- 
dence. 
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1.2 Classical Chaos and Quantum Chaos 

The task of identifying criteria for chaos in quantum mechanics in a manner that extends 
or generahzes the criteria for chaos in classical mechanics meets with several immediate 
difíiculties. Classical chaos is characterized by "extreme sensitivity to initial conditions." 
This property is normally identified with an exponential divergence, on average, in the 
separation between initially nearby classical trajectòries. The rate of this divergence is 
characterized by a set of exponents, called Lyapunov exponents. In the simplest case of 
a mapping, x(n + 1) = F(x(n)), the largest such exponent, \l, may be defined in the 
following way. Let d(n) = Xa(n) — Xb(n) designate the difference between the phase space 
coordinates, and x^, of two nearby classical trajectòries. For small enough |d|, the growth 
of this difference vector is governed by the tangent matrix, 



The eigenvalues of this matrix characterize the local stability of the mapping about the 
trajectory x(n) (6|, g. (If the mapping describes a Hamiltonian system, then these eigen- 
values must sum to zero.) In particular, the Lyapunov exponents are obtained from the 
geomètric mean of these matrices evaluated along some reference trajectory. The largest of 



A more general definition (one that may be applied for practical computation) is provided 
in Chapter 2. 

However, there is no definite concept of "trajectory" in quantum theory. Each quantum 
state, if pure, is specified by a vector in Hilbert space, and not a point in phase space. 
Moreover, the time-dependent Schrodinger equation is linear (more generally, the quantum 
time-development is given by a unitary transformation) , whereas classical chaos requires a 
nonlinear tranformation of the coordinates. For bounded systems, the eigenvalue spectrum 
is discrete, and the quantum dynamics are at most quasi-periodic. 

The significance of these differences may be illustrated by considering some measure of 
proximity that applies to quantum states, and then determining how this measure evolves 
in time. One definition of proximity that may be constructed in Hilbert space is the degree 
of overlap between two quantum states, KV'I^)]. Nearby states are then those that the 




(1.12) 



these exponents may be defined from the double limit |^0| 
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satisfy, KV'I'/')! = 1 — (5 with 6 <^ 1. However, this scalar quantity is invariant under unitary 
time-evolution, 

\{m\m)\ = \{m\ u-\t) u{t) \m)\ = mmmv (i-^) 

Therefore, any two quantum states do not separate at all with time, for any system! 

This attempt at characterizing chaos in quantum mechanics, though unsuccessful, draws 
attention to two of the central qüestions that have motivated research in the field of "quan- 
tum chaos." The first of these pertains to identifying signatures of chaos that arise in 
quantum mechanics: which criteria, if any, distinguish chaotic systems from regular ones 
in quantum theory? The second pertains to quantum-classical correspondence: how does 
classical chaos emerge from the underlying quantum description in the macroscopic limit? 
In particular, how can quantum theory reproduce the exponential divergence exhibited by 
the classical trajectòries! 

Much of the work that attempts to identify manifestations of chaos in quantum theory 
has been obtained from numerical computation of the properties of specific Hamiltonian 
models. These results have suggested a few, possibly universal, signatures of chaos in quan- 
tum systems. Before proceeding to introduce some of these approaches, it is important to 
clarify that there is currently no definitive criterion for chaos in quantum theory. Conse- 
quently, in this thesis, when the term "chaotic" is invoked to describe a quantum system, 
this should be understood as a shorthand expression denoting that the classical counterpart 
to the quantum system exhibits chaotic behaviour. 

Perhaps the most well-known signature of "quantum chaos" arises in the distribution of 
spacings between consecutive energy eigenvalues. Quantum systems with chaotic classical 
counterparts exhibit the same characteristic distributions for these spacings as those of a 
random matrix with the same symmetry properties as the Hamiltonian. These distributions 
exhibit level repulsion relative to their integrable counterparts: small spacings are suppressed 



(leading to "avoided crossings") and most spacings are clustered about an average value |55]. 
This signature exhibited in the eigenvalue spectrum has no obvious classical analogue. 

Although there is no sensitivity to changes in the initial conditions in quantum mechanics 
(given the same unitary evolution for both initial conditions), the presence of classical chaos 
is evident in the quantum state's sensitivity to small perturbations in the Hamiltonian. 
This distinctive feature of the chaotic quantum dynamics has been demonstrated in several 
model systems, although the perturbations have been applied using very different techniques 
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75, Q^, 84, 85, 10 1. The bàsic idea is to evolve a given initial state, |^(0)), according to 
a unitary transformation, |^(í)) = U\^{0)), where U = exp{—iHt/h) for some Hamiltonian 
H, and then examine the overlap, {{ip' {t)\ip{t))\, defined between this final state with another 
perturbed final state, \ip'it)) = C/'|^(0)). Here U' may differ from U according to some small 
perturbation in the Hamiltonian, or by some small coordinate transformation applied at the 



end of the unpertubed evolution [10|. If the classical Hamiltonian is chaotic, then the overlap 



decreases much more rapidly, and to a lower mean asymptotic value, than if the classical 



Hamiltonian is regular [75|. As noted by several authors [p2| , 10 1, this sensitivity criterion 
arises also in the case of a classical density that is evolved according to the Liouville equation. 
As will be demonstrated further below, much can be learned about quantum mechanics, and 
quantum chaos, from the connection between quantum mechanics and the classical Liouville 
equation. 

1.3 Chaos and Ehrenfest Correspondence 

Several authors have addressed the question of quantum-classical correspondence for chaotic 
Systems by comparing the dynamical behaviour of the quantum expectation vàlues with that 



of the classical dynamical variables satisfying Hamilton's equations of motion 48, p 



This approach is clearly motivated by the Ehrenfest correspondence conjecture explained 



in the previous section (1.7). On this view, the quantum dynamics may be identified as 
chaotic provided the quantum expectation vàlues approximately follow a chaotic classical 
trajectory. If nearby classical trajectòries are diverging, on average, exponentially, then the 
quantum expectation vàlues will exhibit a similar exponential divergence, at least while the 
quantum state remains well-localised, that is, in the initial Ehrenfest regime, t < tEhr- 

Accordingly, one of the main interests has been to determine how the correspondence 
time, tEhr, increases with increasing system size. In the literature this limit is often denoted 
using the short-hand expression "?i — > 0", where in this context h always denotes the 
dimensionless ratio {h/J'), and the magnitude of J is given by the quantum numbers or 
some combination of the system parameters. It is also useful to note that {h/J') does not 
actually denote a mathematical limit in quantum theory, but refers to a sequence of quantum 
models with increasing vàlues of the ratio {J/h), but with other parameters adjusted so 
that each quantum model is associated with the same classical model. 

Haake et al. have stuided this time-scale in the special case of a simple non-autonomous 
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spin model which they call the kicked top |54]. The length of the spin is conserved, so 
the quantum number j is a measure of the system dimension jh ~ For initial SU(2) 
coherent states, they obtain an approximate mapping relation for the expectation vàlues 
of the angular momentum components, (Jfc), that consists of the corresponding classical 
mapping modifed by correction terms in powers of (1/j). They are able to relate the growth 
of the first-order corrections to the tangent matrix of the classical mapping. The eigenvalues 
of this matrix govern the local stability of the classical motion. For mappings, the definition 
of the Lyapunov exponent is obtained from the geomètric mean of these matrices evaluated 
along the classical trajectory (see Chapter 2). If the classical motion is regular, then the 
perturbation should grow, at most, as some power of the time. Consequently, for large j, 
significant deviations from the classical map occur on a very late time-scale, 

(1.15) 

where a ~ 1. This scaling is consistent with the example of the baseball considered above. 
On the other hand, if the classical mapping is chaotic, then the first-order quantum correc- 
tions terms should grow, on average, exponentially with time. For chaotic dynamics, the 
quantum mapping therefore deviates from the classical one on the time-scale, 

t^X^hniJ/h). (1.16) 

The above correspondence scaling results for regular ( p..l5 ) and chaotic ( |1.16| ) motion 
are consistent with estimates derived by other authors in several different model systems 



[[L5 , lÉ, 95, 24, 4^, though under a variety of different approximation techniques. In Chapter 
5, I will provide a simpler and much more general argument than the one outlined above, 
which demonstrates that Eqs. ( 1.15| ) and (1.16) should hold (for regular and chaotic systems 
respectively) in the general case of a system of partides interacting only through position 
dependent potentials. 

There are a few features of the above break-time rules that are worth stressing im- 
mediately. First, Eqs. ( |1.15 ) and ( 1.16| ) apply specifically in the case of initial minimum 
uncertainty states, and in this sense they maximize the time-scale of correspondence. Sec- 
ond, these scaling results are guaranteed to apply only when considering the onset of very 
small differences, since the exponential growth of the differences is derived on the basis of 
the local stability of the classical flow. The scaling for large differences, that is, differences 
given by a significant fraction of the system dimension, is left undetermined. Finally, under 
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the restrictive conditions noted just above, these expressions may be taken as an estimate of 
how the Ehrenfest break-time, defined by ( |1.9| ), scales with increasing system size for chaotic 
systems. This follows from the fact that the coordinate centroids of the quantum states are 
associated with the phase space coordinates of a classical trajectory. Consequently, although 
the time-scale ( [L.16 ) is usually cahed simply the "break-time" of quantum-classical corre- 



spondence for a chaotic system, in this thesis I will refer to this result more specificany 
as the Ehrenfest break-time, for reasons that will become clear shortly. As a result of the 
time-scale (|1.16| ), several authors have argued that deterministic chaos can exist only as a 
transient phenomena in quantum mechanics [p2| , [25[| . 

Another distinct feature of quantum chaos was first identified by Fox and coworkers 
in numerical studies of the kicked top and the kicked rotor (which they call the kicked 
pendululm) p7|, |4^. These authors showed that, when the classical dynamics are chaotic, the 
widths of initially localised quantum states grow exponentially with time, until saturation 
at the system size. Specifically, their numerical studies showed that, 

Aq{t) ~ Aq{0) exp X^t, (1.17) 

where Aq'^ = {q^) — {q)"^ is the variance of the quantum state and is a characteristic 
exponent of order the largest Lyapunov exponent. 

A separate, though equally important, contribution of their work was to recognize that 
the exponential spreading of the quantum state mimicks that of a classical distribution 
of trajectòries concentrated around the initial quantum centroids. More precisely, they 
constructed initial classical densities which matched the moments of the Husimi phase space 



functions for the initial coherent states |46|. This close correspondence in the growth rates 



of quantum states and classical densities has been characterized more carefully in recent 



work by Ballentine and McRae |11, 12], who have studied this problem in the Henon-Heiles 
model. The exponential growth of quantum and classical variances is demonstrated also for 
the chaotic states of the coupled spin model examined in this thesis (see Chapter 3). 

This exponential spreading of quantum states is closely related to the earlier estimate 
( |1.16| ) that Ehrenfest correspondence is governed by a \og{J /K) break-time. As noted above, 
the width of the quantum state gives rise to the first correction term in the expansion ( |1.3[ ) , 
and the precision of the correspondence between the quantum expectation vàlues and the 
classical trajectory is limited by the magnitude of this correction term. In particular, for 
an initial minimum-uncertainty state, Aq{tí) ~ ?i/Ap(0), the time at which the width in 
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position space has grown sufficiently large that (^]^) no longer holds is obtained by setting 
the RHS of (1.17) equal to the RHS of ( |1.5D and solving for the Ehrenfest break-time. This 
gives, 

,~A-l„(í^), (L18) 

where L is a characteristic length defined relative to the system size. The magnitude L 
may be interpreted as a threshold indicating a "break" between the quantum and classical 



predictions. It is important to stress that the vahdity of (1.18|) is based on a numerical result 



(1.17), rather than derived, but is more general than the earlier estimate of the break-time 
scahng ( |1.16| ), since the derivation of Eq. ( 1.16| ) apphed strictly in the hmit of vanishingly 
smaU differences. As noted earher, the results of Fox et al. indicate that this exponential 



growth persists until the state width reaches the system size; therefore the scaling rule ( 1.18 ) 
remains vahd even for estimating the onset of very large quantum-classical differences. 

Clearly, for a macroscopic system, the ratio /h is astronomically large. The incredibly 
large magnitude of this ratio is a crucial ingredient in Ehrenfest 's approach to extracting 
Newton's Laws from quantum mechanics However, Zurek and Paz have observed 
that a log(i7//i) time-scale may be short even for systems of macroscopic dimension, i.e., sys- 
tems with J^/h >> 1. Specifically, they have estimated that the time-scale t ~ A^^^ \n{J' /h) 
is comparable to experimentally relevant time-scales in the case of Hyperion, one of the 



moons of Saturn, which is known to exhibit a chaotic tumble |94, Therefore the Ehren- 
fest correspondence principle is at risk of breaking down in the case of chaotic motion. I will 
examine this possibility in detall in Chapter 5, and conclude that such a breakdown actually 
may be subject to experimental confirmation! First, however, in view of the far-reaching 
implications of this result, it is worthwhile to reconsider the assumptions about quantum- 
classical correspondence that are reflected in the Ehrenfest correspondence principle and 
explore more carefully the alternatives to this correspondence conjecture. 



1.4 Quantum Correspondence with Liouville Mechanics 

In the preceding discussion I have drawn attention to some similarities between the prop- 
erties of a time-evolved quantum state and those of a corresponding ensemble of classical 
trajectòries. The general description of the time-evolution of this ensemble is prescribed 
by the classical Liouville equation. This connection may provide further insight into the 
problem of identifying the distinctive features of the quantum dynamics that arise under 
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conditions of chaos in the classical Hamiltonian dynamics. Moreover, in view of the prac- 
tical and conceptual difficulties that challenge the validity of the Ehrenfest correspondence 
conjecture (that is, the view that the predictions of quantum mechanics coincide with New- 
ton's laws for macroscpic objects), it is important to determine whether the correspondence 
between quantum mechanics and classical Liouville mechanics is constrained by similar dif- 
ficulties in the case of chaotic dynamics. 

A useful starting point for investigating the correspondence between these two theories 
starts with a description of the classical Liouville equation. Since one typically does not 
have exact knowledge of the initial coordinates that describe a physical system, Liouville was 
interested in describing how this cloud, or ensemble, of possible initial coordinates evolves 



in time. The time-dependence is given by the following partial differential equation |51|, 

^^ = Mx),ií(x)} (L19) 

where /o(x, í) is a non-negative distribution function (a probability density), x are the phase 
space coordinates, and the symbols {•,•} denote the Poisson bracket. A unique solution 
to the time-dependence requires the specification of an initial boundary condition for the 
density, i.e., some initial state p{q,p, 0). A classical dynamical system is also associated with 
an invariant (Liouville) measure /^(x) on the phase space manifold V. Then, at each point in 
time, one can calculate the ensemble average for an observable j4(x) from the prescription, 

(A(x)), = / dfi{x) /)e(x) A(x). (1.20) 

Jv 

Ballentine et al. have shown that, in the genèric case of Hamiltonians of the form, 
H = /2m+V {q), the time-dependence of the classical ensemble averages exhibit a striking 
similarity to the quantum ones. The Liouville description even includes some of the ostensi- 
bly "quantum" corrections to the classical canonical equations of motion. From (|1.19| ) and 
( 1.201 ) and integrating by parts, one gets, 

d{q)c ^ {p)c 
dt m 

^ = {F{q))c. (1.21) 

As before, it is useful to expand the force function about the position centroid (of the 
classical state), which gives, 

d{q)c ^ {p)c 
dt m 
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dt ' 2 dq^ 



+ ... (1.22) 

q={q)c 



where Ag^ = (g^)c — {q)c· noted above for the time-development of the quantum ex- 
pectation vàlues (1.1), the classical ensmble averages also exibit a first-order "correction" 
term which is proportional to the state variance and which will lead to deviations from the 
canonical equations of motion ( |1.4D . 

Of course, the quantum and classical Liouville theories do not provide identical the- 
oretical descriptions of the same physical system. There are two independent sources of 
discrepency: dynamical differences and differences in the initial conditions. The dynamical 
differences become apparent only when the time-development of the higher-order moments 
(that is, the additional terms in the expansion (1.3)) are calculated explicitly. These calcu- 
lations are carried out most easily in the Wigner-Weyl representation of quantum mechanics 
[p3| . The relevent features of this representation are introduced in Appendix B. 

The time-dependence of the Wigner quasi-distribution pwÍQtP) can be expressed as |71 



dt ^^'^^"^^5^ Uj (2n + l)! V'^+i 5g2n+i 



:i.23) 



where {•, •} denotes the classical Poisson bracket. As explained in Appendix B, although the 
"Moyal terms" contain an explícit proportionality to powers of h, these terms are not guar- 
anteed to become vanishingly small in the limit "7í — > 0" since the momentum derivatives 
will generally produce factors proportional to {q/K)- 

The Wigner quasi-distribution is not a unique quantum phase space distribution, and 
more importantly, may take on negative vàlues. Consequently, it may not be interpreted as 
a probability density and does not have direct experimental significance. The origin of the 
dynamical differences between the quantum and classical theories is clearer if we consider 
differences arising explicitly at the level of observable quantities. As shown in Appendix B, 
for a general time- independent Weyl-ordered operator, A{q,p), the time-development of the 
corresponding quantum expectation value {A{q,p)), is given by, 

d{A^{q,p)) _ , ^fhy- g2n+l^g2n+lff 

dt + (^2j (2n + l)! 5p2-+i 5g2-+i • ^'■^^> 

The Moyal corrections to the dynamics for Liouville averages do not arise explicitly unless 
the operator A{q,p) contains at least a cúbic power of the momentum. Differences arise 
implicitly, however, even for the low-order moments, since the time-development of the 
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low-order moments will eventually become sensitiva to the Moyal corrections through their 
dependence on the time-development of the higher moments. The quantum and Liouville 
time-evolution equations evidently exhibit a much greater degree of similarity than that 
observed between the quantum equations and Hamilton's equations of motion. 

Similary, the Liouvihe picture provides a much more versatile framework for matching 
the details of the initial quantum state. Using a Liouvihe density with non-vanishing sup- 
port, it is possible to match not only the initial quantum expectation vàlues (associated 
with the classical dynamical variables), but also the higher-order quantum moments. The 
"best-case" scenario arises for the coherent states, which are minimum uncertainty states, 
AqAp = h/2. These states are the unique quantum states that are associated with non- 
negative Wigner functions. Consequently all of the quantum moments (given by Hermitian 
Weyl-ordered operators) can be matched exactly by those of a classical density. The appro- 
priate classical density is a Gaussian with support concentrated over an area AçcApc — ^, 
as expected from the uncertainty principle. Perhaps not surprisingly, this matching may 
not be carried out to arbitrary precision for a general quantum state. An intuitive way 
of understanding this is to note that the Wigner quasi-distribution will generally take on 
negative vàlues, and therefore some of the quantum moments will take on numerical vàlues 
that cannot be reproduced by a non-negative Liouville density. 

These formal observations motivate an alternative interpretation of the requirements of 
the correspondence principle. This interpretation consists of the conjecture that the the pre- 
dictions of quantum mechanics should coincide with those of classical Liouville mechanics 
when the predictions of the latter are experimentally vàlid. It will sometimes be conve- 
nient to refer to this conjecture as the Liouville correspondence principle. This is a weaker 
conjecture than the Ehrenfest correspondence principle, which was stated in the discussion 
preceding the heurístic formula (1.7). 

Ballentine and coworkers have studied the details of this quantum-Liouville correspon- 
dence for a few non-integrable systems. Ballentine, Yang, and Zibin compared quantum 
expectation vàlues and classical ensemble averages for the low-order moments, e.g. the dy- 
namical variables and their variances, primarily at early times, that is, in the initial Ehrenfest 
regime of a few low-dimensional model systems. For fixed system size, the correspondence 
with Liouville mechanics was shown to be much more accurate, and to last much longer, 
than the correspondence with a single classical trajectory. They demonstrated a remarkable 
degree of improved correspondence under conditions of classical chaos. 
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Ballentine and McRae |11| and Ballentine |jï^ examined the degree of Liouville corre- 
spondence for the Henon-Heiles model. They demonstrated that both the quantum and 
classical variances grow at an approximately exponential rate, a result obtained previously 
by Fox and coworkers for different model systems They also characterised the dy- 

namical behaviour of the quantum-LiouvilIe differences for the dynamical variables and their 
variances. They determined that, for states initiated from a chaotic region, the quantum- 
Liouville differences grow exponentially, in the initial Ehrefenst regime, with a characteristic 
exponent that is larger than the largest Lyapunov exponent. 

The duration of the correspondence between quantum mechanics and Liouville mechanics 
has been considered by Zurek and Paz [ p6| , p^ ] in the "optimal" case of minimum uncertainty 
states. Using very general arguments about the nature of the chaotic classical dynamics and 
the growth of the Moyal corrections in Eq. ( |1.23| ), they have estimated that the quantum- 
Liouville differences should be governed by a break time that scales as, 

~ -3-ln(L Ap(0)/n) (1.25) 
Al 



where L and Ap(0) are defined as in Eqs. ( |1.5| ) and (1.18). Casati and Chirikov [^] have 
arg ued that ([L25| ) is not a distinct result, but is essentially the same as the (Ehrenfest) break- 
time, given by ( 1.16] ), that was first developed in Refs. [Ï5|, 24]. Although the estimates ( |1.25 ) 
and ( |1.16|) bear a formal similarity, this Introduction makes clear that they are based on 
distinct correspondence criteria. (The conditions for Liouville and Ehrenfest correspondence 
are often not clearly distinguished in the literature.) 

However, the formal similarity between tb and ÍEhr suggests that the Liouville corre- 
spondence principle may also be at risk of breaking down in the case of chaotic macroscopic 
motion, for systems such as Hyperion, since the break-time scaling exhibited by Eq. ( 1.251 ) 
suggests the possibility of observable quantum-classical differences on an experimentally 
accessible time-scale p7| . While Zurek has raised the possibility of a breakdown of the cor- 
respondence principle |^9| in this context, he has also suggested that this breakdown is not 
actually observed because of the "decoherence" effects that result from weak coupling to 
the ubiquitous environment. 

Habib et al. have considered the effects of decoherence on the degree of Liouville cor- 
respondence at the level of expectation vàlues for low-order moments and also by comparing 
Wigner quasi-distributions directly with the classical density in study of the Henon-Heiles 
model. They demonstrated that the degree of correspondence (for fixed h/J') improved 
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when the effects of a stochastic enviroment (or "decoherence" ) are taken into account. 
The positive effects of decol·ierence on the degree of quantum-Liouville correspondence 
have been demonstrated in several recent theoretical studies of classicahy chaotic models 



p2[ p8| , 18, In the case of the kicked rotor model, a suppression of quantum localisa- 
tion due to environmental noise has been demonstrated experimentally in a simulation of 
the model using trapped cesium àtoms subjected to a pulsed làser beam Q. The "envi- 
ronment" in this case consists of spontaneous emission due to the "vacuum fluctuations." 



However, Farini et al. |37] have shown that the magnitude of the Moyal terms in Eq. ( 1.23 ) 
can remain small even for an isolated system, and therefore the effects of decoherence may 
not be necessary to recover classical properties. 



The scaling of the Liouville break-time estimate ( 1.25 ) has been examined by Roncaglia 
et al. [79| in numerical studies of the anomolous dijfusion regime of the kicked rotor. These 
authors provided some numerical evidence in support of the scaling with increasing system 
size expressed in Eq. ( |1.25| ) by considering the differences between the quantum expectation 
vàlues and Liouville ensemble averages for low-order moments, such as (p^). However, it 
is important to emphasize that the log{J^/h) scaling for the Liouville break-time does not 
hold in the general case. In particular, the quantum-Liouville break-time that arises in the 
localisation regime of the kicked rotor is believed to grow as a small power of the ratio {J /fi) 
[^5| , which is more than adequate for correspondence in the macroscopic limit. 

These considerations have been outlined in order to motivate the following qüestions: 
(1) Under what conditions does the Liouville break-time estimate ( |1.25 ) hold? (2) Although 
environmental noise may improve the degree of quantum-Liouville correspondence, is deco- 
herence necessary for quantum-Liouville correspondence in the macroscopic limit? (3) Can 
decoherence help restore Ehrenfest correspondence? 

The main focus of the work presented in this thesis is to characterize the properties of 
Liouville and Ehrenfest correspondence, particularly under conditions of classical chaos, for 
an isolated few degree-of-freedom system. The central challenge of this analysis consists 
of extrapolating the predictions of quantum theory (and therefore the degree of correspon- 
dence) into the macroscopic regime, where direct calculation is of course not feasible in the 
genèric case of a non-integrable system. 

This thesis is organised as follows. In Chapter 2 I will describe the quantum and classical 
versions of the model examined in this thesis. Since the model is novel I will examine the 
behaviours of the chaotic classical dynamics in some detail. I will also describe the initial 
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quantum states, which are SU(2) coherent states, and then define a corresponding classical 
density on the 2-sphere which is a good analog for these states. I show in Appendix A that 
a perfect match is impossible: no distribution on can reproduce the moments of the 
SU(2) coherent states exactly. At the end of Chapter 2 I will describe some of the numerical 
techniques that were required for the quantum and classical Liouville descriptions of the 
dynamics. 

In Chapter 3 I will examine the characteristics of quantum-Liouville correspondence at 
the level of the low-order moments. First I will demonstrate that a close correspondence 
with Liouville mechanics persists well after the Ehrenfest correspondence has broken down. 
The widths of both the quantum states and the Liouville densities grow exponentially until 
saturation at the system size. I will then confirm that the quantum-Liouville differences 
exhibit an initial phase of exponential growth for chaotic states, but also demonstrate that 
this growth terminates well before the differences reach the system size. I will also show 
that the Liouville break-time scales according to the estimate Eq. ( 1.25| ), but only for the 
onset of very small quantum-Liouville differences. By "very small" I mean differences that 
remain a factor {h/J') smaller than the system size. As a result, I argue that decoherence 
is not necessary for the emergence of these classical properties in the macroscopic limit. 

This demonstrated correspondence at the level of the low-order moments still leaves 
room for significant quantum-Liouville deviations. In Chapter 4 I will examine the degree 
of correspondence between the quantum and classical probability distributions, emphasiz- 
ing the detalls of this correspondence after the states have spread to the system size. In 
particular, I will show that the initially localised quantum states approach the classical 
equilibrium on the same relaxation time-scale as the corresponding Liouville density. When 
the classical Hamiltonian exhibits widespread chaos, the quantum probability distributions 
approach microcanonical equilibrium configurations. I also demonstrate that the quantum 
fluctuations away from the classical equilibrium become vanishingly small in the macro- 
scopic limit. These results reinforce the conclusion that quantum-Liouville differences for 
classical observables become increasingly difhcult to detect as the system size increases. 

In Chapter 5 I will examine the properties of Ehrenfest correspondence. I will demon- 
strate that the Ehrenfest differences, in contrast with the quantum-Liouville differences, 
actually grow to a macroscopic size on the short time-scale ( p..l8D in the case of chaotic clas- 
sical dynamics. Moreover, decoherence is unable to restore the Ehrenfest correspondence 
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for these observables. I will conclude that Ehrenfest correspondence may be subject to ex- 
perimental refutation for some chaotic macroscopic bodies. At the end of Chapter 5 I will 
address the implications of this result for the correspondence principle and the interpretation 
of the quantum state. 



Chapter 2 



Description of the Model 



This chapter is taken from Emerson and Ballentine [35|. 



2.1 The Model 



We consider the quantum and classical dynamics generated by a non-integrable model of 
two interacting spins, 



where S = {Sx, Sy, Sz) and L = {L^, Ly, Lz). The first two terms in (|2.1|) correspond to 
simple rotation of both spins about the z-axis. The sum over coupling terms describes an 
infinite sequence of 5-function interactions at times t = n for integer n. Each interaction term 
corresponds to an impulsive rotation of each spin about the j;-axis by an angle proportional 
to the x-component of the other spin. 

2.1.1 The Quantum Dynamics 

To obtain the quantum dynamics we intèrpret the Cartesian components of the spins as 
operators satisfying the usual angular momentum commutation relations, 



oo 
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In the above we have set fi = 1 and introduced the total angular momentum vector J = S+L. 

The Hamiltonian ( |2.1| ) possesses kinematic constants of the motion, [S^,-fí] = and 
[L^,ií] = 0, and the total state vector \tp) can be represented in a finite Hilbert space of 
dimension N = {2s + 1) x [21 + 1). This space is spanned by the orthonormal vectors 
|s, l, nisjUii) = \s, nis) |/, mi) with rUs £ {s, s — 1, . . . , —s} and mi £ {1,1 — 1, ... , —l}. 
These are the joint eigenvectors of the four spin operators 

S'^\s,l,ms,mi) = s{s + l)\s,l,ms,mi), 

Sz\s,l, 171^,1711) = ■ms\s,l,ms,mi), (2.2) 

'L'^\s,l,ms,mi) = l{l + l)\s,l,ms,mi), 

Lz\s,l,ms,mi) = mi\s,l,ms,mi). 

The periòdic sequence of interactions introduced by the (5-function produces a quantum 
mapping. The time-evolution for a single iteration, from just before a kick to just before 
the next, is produced by the unitary transformation, 

\^ÍJ{n + l))=F\^l;{n)), (2.3) 

where F is the single-step Floquet operator, 

F = exp [-ia{Sz + L;,)] exp [-icSxL^] . (2.4) 

Since a is a rotation its range is 27r radians. The quantum dynamics are thus specified by 
two parameters, a and c, and two quantum numbers, s and /. 

An explícit representation of the single-step Floquet operator can be obtained in the basis 
( |2.2D by first re-expressing the interaction operator in ( |2.4| ) in terms of rotation operators, 

exp [-icSx Lx] = (6», 0) R^^^ {9, cj))] exp [-icS^ L,] 

X [i? W(e,</.)®iï«(e, </.)]-!, (2.5) 

using polar angle 9 = Tr/2 and azimuthal angle (p = 0. Then non-zero off-diagonal terms 
arise only in the expressions for the rotation matrices, which take the form, 

{j,m'\R^^\9,^)\j,m) =exp{-im'^)d'£,J9). (2.6) 

The matrix elements, 

€',m(G) = eM-iOJy)\j,m) (2.7) 
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are given explicitly by Wigner's formula [31|. 

We are interested in studying the different time-domain characteristics of quantum ob- 
servables when the corresponding classical system exhibits either regular or chaotic dynam- 
ics. In order to compare quantum systems with different quantum numbers it is convenient 
to normalize subsystem observables by the subsystem magnitude ^/{Lt^ = \JI{1 + 1). We 
denote such normalized observables with a tilde, where 



(4(n)) 



{21 



and the normalized variance at time n is defined as, 

We are also interested in evaluating the properties of the quantum probability distri- 
butions. The probability distribution corresponding to the observable is given by the 
trace, 



P^(m/) = Tr p^^\n)\l,mi){l,mi\ = {l,mi\p^^^{n)\l,mi), 



(2.10) 



where , rrisl ] is the reduced state operator for the spin 

L at time n and Tr^**) denotes a trace over the factor space corresponding to the spin S. 



2.1.2 Classical Map 

For the Hamiltonian (^]^) the corresponding classical equations of motion are obtained by 
interpreting the angular momentum components as dynamical variables satisfying, 



{Si, Sj} 

{L^,LJ} 

{Ji,jj} 



^ijk-^k 
^ijkJki 



with {•, •} denoting the Poisson bracket. The periòdic (5-function in the coupling term can 
be used to define surfaces at í = n, for integer n, on which the time-evolution reduces to a 
stroboscopic mapping, 



■n+l 



òn+l 



S'"cos(a)- S'^cos(7rL") - 5'"sin(7rL") sin(a), 
Sy cos(7rL") - 5" sin(7rZ") cos(a) + 5" sin(a). 
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'z 
n+1 

X 

n+1 

y 



L 



L 



4"cos(7rLS)+5,"sin(7rL,"), 
L"cos(a)— Ly cos{-ySl^) — sm{'yS^) sm(a), 
cos(75") - L"sin(7S'") cos(a) + Z"sin(a), 
L^cos(75^)+L^sin(75^), 



(2.11) 



where L = L/|L| , S = S/|S| and we have introduced the parameters 7 = c|S| and r = 
|L|/|S|. The mapping equations ( 2.11 ) describe the time-evolution of (|2.l| ) from just before 
one kick to just before the next. 

Since the magnitudes of both spins are conserved, {S^,H} = {L^^H} = 0, the motion 
is actuaUy confined to the four-dimensional manifold V = S'^ x S^, which corresponds 



to the surfaces of two spheres. This is manifest when the mapping (2.11) is expressed 
in terms of the four canonical coordinates x = {S^, (ps, L^, (pi), where (ps = ^Si.n(Sy/ Sx) and 
(pl = tan^Ly/Lx). We wih refer to the mapping (|2.11| ) in canonical form using the shorthand 
notation x""*'^ = F(x"'). It is also useful to introduce a complete set of spherical coordinates 
e = {es,(ps,9i,(pi) where 0, = cos-i(5./|S|) and 61 = cos-i(L,/|L|). 

The mapping ( |2.11| ) on the reduced surface V enjoys a rather large parameter space; 
the dynamics are determined from three independent dimensionless parameters: a £ [0, 2tt), 
7 G (—00, 00), and r > 1. The first of these, a, controls the angle of free-field rotation about 
the z-axis. The parameter 7 = c|S| is a dimensionless coupHng strength and r = |L|/|S| 
corresponds to the relative magnitude of the two spins. 

We are particularly interested in the effect of increasing the couphng strength 7 for 



different fixed vàlues of r. In Fig. 2.1 we plot the dependence of the classical behaviour 



on these two parameters for the case a = 5, which produces typical results. The data in 
this figure were generated by randomly samphng initial conditions on V, using the invariant 
measure, 

dii{x) = dSzd(psdLzd(pi, (2-12) 

and then calculating the largest Lyapunov exponent associated with each trajectory. Open 
circles correspond to regimes where at least 99% of the initial conditions were found to 
exhibit regular behaviour and crosses correspond to regimes where at least 99% of these 
randomly sampled initial conditions were found to exhibit chaotic behaviour. Circles with 
crosses through them (the superposition of both symbols) correspond to regimes with a 
mixed phase space. For the case a = 5 and with r held constant, the scaled coupling 
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Figure 2.1: Behaviour of the classical mapping for different vàlues of r = |L|/|S| and 
7 = c|S| with a = 5. Circles correspond to parameter vàlues for which at least 99% of the 
surface area V produces regular dynamics and crosses correspond to parameter vàlues for 
which the dynamics are at least 99% chaotic. Superpositions of circles and crosses correspond 
to parameter vàlues which produce a mixed phase space. We investigate quantum-classical 
correspondence for the parameter vàlues 7 = 1.215 (mixed regime) and 7 = 2.835 (global 
chaos), with r = 1.1, which are indicated by filled circles. 



strength 7 plays the role of a perturbation parameter: the classical behaviour varies from 
regular, to mixed, to predominantly chaotic as I7I is increased from zero. 

The fixed points of the classical map ( 2.11| ) provide useful information about the pa- 
rameter dependence of the classical behaviour and, more importantly, in the case of mixed 
regimes, help locate the zones of regular behaviour in the 4-dimensional phase space. We 
find it sufíicient to consider only the four trivial (parameter-independent) fixed points 
which lie at the poles along the z-axis: two of these points correspond to parallel spins, 
{Sz,Lz) = =b(|S|,|L|), and the remaining two points correspond to anti-parallel spins, 
(5„L,) = (±|S|,T|L|). 

The stability around these fixed points can be determined from the eigenvalues of the 
tangent map matrix, M = dF/dx, where all derivatives are evaluated at the fixed point of 



interest. (It is easiest to derive M using the six non-canonical mapping equations ( 2.11 ) 
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since the tangent map for the canonical mapping equations exhibits a coordinate system 
singular ity at these fixed points.) The eigenvalues corresponding to the four trivial fixed 
points are obtained from the characteristic equation, 



[r - 2^ cosa + 1]^ ± CTrsm^a = 0, (2.13) 

with the minus (plus) sign corresponding to the parallel (anti-parallel) cases and I have 
suppressed the trivial factor (1 — ^)^ which arises since the six equations ( 2.11 ) are not 
independent. For the parallel fixed points the four eigenvalues are 



^12 = cos a ± -\l r^"^ svo? a + -\j ±4 cos a\l ^"^r siv? a — [svo? a)(4 — 7^r), 



1 / . ■ . 1 



Ç|^4 = COS a ib — Y sin^ a — — Y ib4cos ay 7^r sin^ a — (sin^ a)(4 — 7^r), (2.14) 

and the eigenvalues for the anti-parallel cases, Ç^^, are obtained from ( |2.14| ) through the 
substitution r — > — r. A fixed point becomes unstable if and only if |^| > 1 for at least one 
of the four eigenvalues. 

Mixed Phase Space: 7 = 1.215 

We are particularly interested in the behaviour of this model when the two spins are com- 
parable in magnitude. Choosing the value r = 1.1 (with a = 5 as before), we determined 
by numerical evaluation that the anti-parallel fixed points are unstable for I7I > 0. In the 
case of the parallel fixed points, all four eigenvalues remain on the unit circle, = 1, for 
I7I < 1.42. This stability condition guarantees the presence of regular islands about the 



parallel fixed points ||6q| . In Fig. 2.2 we plot the trajectory corresponding to the parameters 
a = 5, r = 1.1, 7 = 1.215 and with initial condition 9{tí) = (5°, 5°, 5°, 5°) which locates the 
trajectory near a stable fixed point of a mixed phase space (see Fig. |2.1| .) This trajectory 
clearly exhibits a periòdic pattern which we have confirmed to be regular by computing the 
associated Lyapunov exponent (A^ = 0). In contrast, the trajectory plotted in Fig. O is 
launched with the same parameters but with initial condition 0(0) = (20°, 40°, 160°, 130°), 
which is close to one of the unstable anti-parallel fixed points. This trajectory explores 
a much larger portion of the surface of the two spheres in a seemingly random manner. 
As expected, a computation of the largest associated Lyapunov exponent yields a positive 
number (Al = 0.04). 
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Figure 2.2: Stroboscopic trajectòries on the unit sphere launched from a regular zone of the 
mixed regime with 7 = 1.215, r = 1.1, a = 5 and 9{0) = (5°, 5°, 5°, 5°). 

Global Chaos: 7 = 2.835 

If we increase the coupling strength to the value 7 = 2.835, with a = 5 and r = 1.1 as 
before, then all four trivial fixed points become unstable. By randomly sampling V with 
3 X 10"^ initial conditions we find that less than 0.1% of the kinematically accessible surface 
V is covered with regular islands (see Fig. |2.1| ) . This set of parameters produces a connected 
chaotic zone with largest Lyapunov exponent = 0.45. We will refer to this type of regime 
as one of "global chaos" although the reader should note that our usage of this expression 
differs slightly from that in Ref. [p^]. 

The Límit r > 1 

Another interesting limit of our model arises when one of the spins is much larger than 
the other, r 1. We expect that in this limit the larger spin (L) will act as a source of 
essentially external "driving" for the smaller spin (S). Referring to the coupling terms in 
the mapping ( |2.11| ), the "driving" strength, or perturbation upon S from L, is determined 
from the product jr = c|L|, which can be quite large, whereas the "back-reaction" strength, 
or perturbation upon L from S, is governed only by the scaled coupling strength 7 = c|S|, 
which can be quite small. It is interesting to examine whether a dynamical regime exists 
where the larger system might approach regular behaviour while the smaller 'driven' system 
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Figure 2.3: Same parameters as Fig. but the trajectory is launched from a chaotic zone 
of the mixed regime with initial condition (9(0) = (20°, 40°, 160°, 130°). 

is still subject to chaotic motion. 

In Fig. we plot a chaotic trajectory for r = 100 with initial condition 6{0) = 
(27°, 27°, 27°, 27°) which is located in a chaotic zone (A^ = 0.026) of a mixed phase space 
(with a = 5 and 7 = 0.06). Although the small spin wanders chaotically over a large portion 
of its kinematicany accessible shell S'^, the motion of the large spin remains confined to a 
'narrow' band. Although the band is narrow relative to the large spin's length, it is not 
small relative to the smaller spin's length. The trajectòries are both plotted on the unit 
sphere, so the effective area explored by the large spin (relative to the effective area covered 
by the small spin) scales in proportion to r^. 



2.1.3 The Liouville Dynamics 



We are interested in comparing the quantum dynamics generated by (2.3) with the corre- 
sponding Liouville dynamics of a classical distribution. The time-evolution of a classical 
phase space distribution is generated by the partial differential equation, 

5/0c(x,í) 



dt 



-{Pc,H}, 



(2.15) 



where H stands for the Hamiltonian (|2.1|) and x = (Sz, (ps,Lz,(pi) 
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Figure 2.4: A chaotic trajectory for mixed regime parameters 7 = 0.06, r = 100, and a = 5 
with 9(0) = (27°, 27°, 27", 27°). The motion of the larger spin appears to remain confined 
to a narrow band on the surface of the sphere. 

The solution to (^l5|) can be expressed in the compact form, 

Pe(x, t)= í d/i(y) 5{x - x(í, y)) p,iy, 0), (2.16) 
Jv 

with measure dfi{y) given by ( ^.12]) . Each time-dependent function x(t,y) G P is a solution 
of the equations of motion for ( |2.1| ) with initial condition y G "P. The solution ( 2. 16] ) 



simply expresses that Liouville's equation (2.15) describes the dynamics of a classical density 



/9c(x, t) of points evolving in phase space under the Hamiltonian flow. We exploit this fact to 
numerically solve ( |2.15| ) by randomly generating initial conditions consistent with an initial 
phase space distribution pc(x, 0) and then time-evolving each of these initial conditions 
using the equations of motion ( p.llj) . We then calculate the ensemble averages of dynamical 
variables, 

{L,{n))c= í dM(x)-^pe(x,n). (2.17) 
by summing over this distribution of trajectòries at each time step. 

2.1.4 Correspondence Between Quantum and Classical Models 

For a quantum system specified by the four numbers {a, c, s, l}, the corresponding classical 
parameters {a, 7, r} may be determined by first defining the classical magnitudes in terms 
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of the quantum magnitudes, 



S|c = 



L|c = 



(2.18) 



where the quantities on the left hand side are the lengths of the elassical spins and those on 
the right are the quantum numbers. If wc set the Hamiltonian coeficients a and c numeri- 
caUy equal for both models, then the remaining two dimensionless elassical parameters are 
determined, 



We are interested in extrapolating the behaviour of the quantum dynamics in the limit 
s oo and / oo. This is accomplished by studying sequcnces of quantum models with 
increasing s and / choscn such that r and 7 arc held fixed. Sincc s and / arc rcstrictcd to 
integer (or half-intcger) vàlues, the corrcsponding elassical r will actually vary slightly for 
each member of this sequence, although 7 can be matched cxactly by varying the quantum 
parameter c. In the limit s ^ 00 and / 00 this variation becomes increasingly small 
since r = + l)/s(s + 1) l/s. For convenience, the elassical r corrcsponding to each 
member of the sequence of qTiantum models is identificd by its value in this limit. Wc have 
examined the effect of the small variations in the value of r on the elassical behaviour and 
found the variation to be negligible. 

2.2 Initial States 

2.2.1 Initial Quantum State 

We consider initial quantum states which are pure and separable. 



r 



s(s + l) 



7 




(2.19) 



\m) = \Mo))^mo)). 



(2.20) 



For the initial state of each subsystem we use one of the directed angular momentum states, 

\e,<p) = R^^\e,<p)\j,j), (2.21) 
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which correspond to states of maximum polarization in the direction {0,(p). It has the 
properties: 

{e,^\j,\e,<p) = jcos9 

{9,cl)\J,±iJy\0,cl)) = je^'^sinO, (2.22) 

where j in this section refers to either / or s. 

The states ( 2.21| ) are the SU(2) coherent states, which, Uke their counterparts in the 



EucHdean phase space, are minimum uncertainty states ||T^; the normahzed variance of the 
quadratic operator, 

^ j2 ^ {9,<I>\J'\0,<I>)-{0A\J\0A)' ^ 1 .2 231 

is minimised for given j and vanishes in the hmit j oo. The coherent states polarized 
along the z-axis, \j,j) and |j, —j), also saturate the inequahty of the uncertainty relation, 

{J!){j'y) > (2.24) 
although this inequahty is not saturated for coherent states polarized along other axes. 

2.2.2 Initial Classical State and Correspondence in the Macroscopic Limit 

We compare the quantum dynamics with that of a classical Liouville density which is chosen 
to match the initial probability distributions of the quantum coherent state. For quantum 
systems with a Euclidean phase space it is always possible to construct a classical density 
with marginal probability distributions that match exactly the corresponding moments of 
the quantum coherent state. This follows from the fact that the marginal distributions for 
a coherent state are positive definite Gaussians, and therefore all of the moments can be 
matched exactly by choosing a Gaussian classical density. For the SU(2) coherent state, 
however, we show in Appendix A that no classical density has marginal distributions that 
can reproduce even the low order moments of the quantum probability distributions (except 
in the limit of infinite j). Thus from the outset it is clear that any choice of an initial 
classical state will exhibit residual discrepancy in matching some of the initial quantum 
moments. 

We have examined the initial state and dynamical quantum-classical correspondence 
using several different classical distributions. These included the vector model distribution 
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described in Appendix A and the Gaussian distribution used by Fox and Elston in corre- 
spondence studies of the kicked top |46]. For a state polarized along the z-axis we chose the 
density, 



Pc 



sm 



C exp 
C exp 



2sin2(|)) 

{1-Jz 



sin 6d6d(j) 
dJ^dò, 



(2.25) 



with C = [2'K(t'^ (1 — exp(— 2(7 2))] ^ ^ instead of those previously considered, because it is 
periòdic under 27r rotation. An initial state directed along {6o-,(j)o) is then produced by a 



rigid body rotation of ( 2.25 ) by an angle 9(, about the y-axis followed by rotation with angle 
4>o about the z-axis. 

The variance a"^ and the magnitude |J|c are free parameters of the classical distribution 
that should be chosen to fit the quantum probabilities as well as possible. It is shown in 
Appendix A that no classical density has marginal distributions which can match all of the 
quantum moments, so we concentrate only on matching the lowest order moments. Since the 
magnitude of the spin is a kinematic constant both classically and quantum mechanically, 
we choose the squared length of the classical spin to have the correct quantum value, 



{jl)c + {Jl)c + {Jl)c=j{j + l)- 



(2.26) 



For a state polarized along the z-axis, we have {Jx) = {Jy) = and (Jy) = (J^) for 
both distributions as a consequence of the axial symmetry. Furthermore, as a consequence 



of ( 2.26 ), we will automatically satisfy the condition, 

2(4)c + (jf>c = j(j + l). 



(2.27) 



Therefore we only need to consider the classical moments, 

{J,)c = |J| G{a^) 
(4), = |J| V G{a^), 

calculated from the density ( p.25 ) in terms of the remaining free parameter, 172^ where, 

"1 + exp(-2(7-2 



2.28) 
2.29) 



1 — exp(— 2(7 



(7 



(2.30) 
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We would like to match both of these classical moments with the corresponding quantum 
vàlues, 

(Jz) = j, (2.31) 
(4')=iA (2.32) 

calculated for the coherent state ( |2.21| ). However, no choice of a"^ will satisfy both con- 
straints. 

If we choose a"^ to satisfy (2.31) exactly then we would obtain, 



If we choose cj^ to satisfy ( 2.32| ) exactly then we would obtain, 



a^ = ^ + ^ + 0(r^). (2.34) 

(These expansions are most easily derived from the approximation G{a'^) ~ 1 — o"^, which 
has an exponentially small error for large j.) 

We have chosen to compromise between these vàlues by fixing o"^ so that the ratio 
{Jz)c/ {Jx)c has the correct quantum value. This leads to the choice, 

These unavoidable initial differences between the classical and quantum moments will 
vanish in the "classical" limit. To see this explicitly it is convenient to introduce a measure 
of the quantum-classical differences, 

SUn) = \{Un))-{Un)),\, (2.36) 



defined at time n. For an initial state polarised in direction (6, (p), the choice ( 2.35 ) produces 
the initial difference, 

SU0) = '-^ + O{r'), (2.37) 

which vanishes as j — > oo. 
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2.3 Numerical Methods 

The time-periodic spin Hamiltonian ( |2.1| ) is convenient for numerical study because the 
time-dependence reduces to a simple mapping and the quantum state vector is confined to 
a finite dimensional Hilbert space. Consequently we can solve the exact time-evolution 
equations (^|^) numerically without introducing any artificial truncation of the Hilbert 
space. The principal source of numerical inaccuracy arises from the numerical evaluation 
of the matrix elements of the rotation operator {j,m'\R{9, (/))\j,m) = eyi-p{—i(l)m')é'^,^{9). 
The rotation operator is required both for the calculation of the initial quantum coherent 
state, \0^(l)) = R{9 , , m = j), and for the evaluation of the unitary Floquet opera- 
tor. In order to maximise the precision of our results we calculated the matrix elements 
^m'm(^) ~ 0' 6xp(— i^Jj/)|j, m) using the recursion algorithm of Ref. |21] and then tested 
the accuracy of our results by introducing controlled numerical errors. For small quantum 
numbers (j < 50) we are able to confirm the correctness of our coded algorithm by com- 
paring these results with those obtained by direct evaluation of Wigner's formula for the 
matrix elements d^^,^{9). 

The time evolution of the Liouville density was simulated by numerically evaluating 
between 10^ and 10^ classical trajectòries with randomly selected initial conditions weighted 
according to the initial distribution ( p.25 ). Such a large number of trajectòries was required 



in order to keep Monte Cario errors small enough to resolve the initial normalized quantum- 
classical differences, which scale as l/8j^, over the range of j vàlues we have examined. 

We identified initial conditions of the classical map as chaotic by numerically calculating 
the largest Lyapunov exponent, Al, using the formula, 

1 ^ 

^L = j^Y.^r,d[n) (2.38) 

n=l 

where d{n) = |5xj(re)| , with d{Çi) = 1. The differential (5x(n) is a diff'erence vector 
between adjacent trajectòries and thus evolves under the action of the tangent map, 

5x(n + l) = M·òx(n), (2.39) 



where M is evaluated along some fiducial trajectory 

Since we are interested in studying quantum states, and corresponding classical distri- 
butions which have non-vanishing support on the sphere, it is also important to get an idea 
of the size of these regular and chaotic zones. By comparing the size of a given regular or 
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chaotic zone to the variance of an initial state located within it, we can determine whether 
most of the state is contained within this zone. However, we can not perform this com- 
parison by direct visual inspection since the relevant phase space is 4-dimensional. One 
strategy which we used to overcome this difhculty was to calculate the Lyapunov exponent 
for a large number of randomly sampled initial conditions and then project only those points 
which are regular (or chaotic) onto the plane spanned by Sz = cos Og and = cos 6i. If the 
variance of the initial quantum state is located within, and several times smaller than, the 
dimensions of a zone devoid of any of these points, then the state in question can be safely 
identified as chaotic (or regular). 



Chapter 3 

Correspondence for Low Order 
Moments 



This chapter is taken from Emerson and Ballentine [35| 



3.1 Introduction 

There is considerable interest in the interface between quantum and classical mechanics and 
the conditions that lead to the emergence of classical behaviour. In order to characterize 
these conditions, it is important to differentiate two dynamical regimes of quantum-classical 
correspondence [0: 

(i) Ehrenfest correspondence, in which the centroid of the wave packet approximately fol- 
lows a classical trajectory; and 

(ii) Liouville correspondence, in which the quantum probability distributions are in approx- 
imate agreement with those of an appropriately constructed classical ensemble satisfying 
Liouville's equation. 

Regime (i) is relevant only when the width of the quantum state is small compared to 
the dimensions of the system; if the initial state is not narrow, this regime may be absent. 
Regime (ii), which generally includes (i), applies to a much broader class of states, and this 
regime of correspondence may persist well after the Ehrenfest correspondence has broken 
down. The distinction between regimes (i) and (ii) has not always been made clear in the 
literature, though the conditions that delimit these two regimes, and in particular their 
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scaling with system parameters, may be quite different. 

The theoretical study of quantum chaos has raised the question of whether the quantum- 
classical break occurs differently in chaotic states, in states of regular motion, and in mixed 
phase-space systems. This is well understood only in the case of regime (i). There it is 



weh-known |15, 54, ^ that the time for a minimum-uncertainty wave packet to expand 
beyond the Ehrenfest regime scales as log{J^/fi) for chaotic states, and as a power oi J /Ti 
for regular states, where J denotes a characteristic system action. 

The breakdown of quantum-classical correspondence, in the case of regime (ii), is less 
well understood, though it has been argued that this regime may also be delimited by 
a \og{J /h) break-time in classically chaotic states |5^. Some numerical evidence in 
support of this conjecture has been reported in a study of the kicked rotor in the anomolous 
diffusion regime |^^. (On the other hand, in the regime of quantum localization, the break- 
time for the kicked rotor seems to scale as {J /K)^ f^^-) Since the log(^/?i) time scale is 
rather short, it has been suggested that certain macroscopic objects would be predicted to 
exhibit non-classical behaviour on observable time scales |97, These results highlight 
the importance of investigating the characteristics of quantum-classical correspondence in 
more detall. 

In this paper we study the classical and quantum dynamics of two interacting spins. This 
model is convenient because the Hilbert space of the quantum system is finite-dimensional, 
and hence tractable for computations. Spin models have been useful in the past for exploring 
classical and quantum chaos ^ ^, |^, |7^ and our model belongs to a class of spin 



models which show promise of experimental realization in the near future |6£]. The classical 
limit is approached by taking the magnitude of both spins to be very large relative to fi, 
while keeping their ratio fixed. For our model a characteristic system action is given by 
J ~ where / is a quantum number, and the classical limit is simply the limit of large 
quantum numbers, i.e. the limit / — > oo. 

In the case of the chaotic dynamics for our model, we first show that the widths of 
both the quantum and classical states grow exponentially at a rate given approximately 
by the largest Lyapunov exponent (until saturation at the system dimension). We then 
show that the initially small quantum-classical differences also grow at an exponential rate, 
with an exponent Açc that is independent of the quantum numbers and at least twice as 
large as the largest Lyapunov exponent. We demonstrate how this exponential growth 
of differences leads to a log break-time rule, th ~ \~^hi{lp/fi), delimiting the regime of 
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Liouville correspondence. The factor p, measured in units of h, is some preset tolerance that 
defines a break between the quantum and classical expectation vàlues, but does not scale 
with the system dimension. However, we also show that this logarithmic rule holds only if 
the tolerance p for quantum-classical differences is chosen extremely small, in particular p < 
0{h). For larger vàlues of the tolerance, the break-time does not occur on this log time-scale 
and may not occur until the recurrence time. In this sense, log break-time rules describing 
Liouville correspondence are not robust. These results demonstrate that, for chaotic states in 
the classical limit, quantum observables are described approximately by Liouville ensemble 
averages well beyond the Ehrenfest time-scale, after which both quantum and classical states 
relax towards equilibrium configurations. This demonstration of correspondence is obtained 
for a few degree-of-freedom quantum system of coupled spins that is described by a pure 
state and subject only to unitary evolution. 



3.2 Characteristics of the Quantum and Liouville Dynamics 
3.2.1 Mixed Phase Space 



We consider the time-development of initial quantum coherent states ( 2.21j ) evolved ac- 



cording to the mapping (|2.3| ) using quantum numbers s = 140 and / = 154 and associ- 
ated classical parameters 7 = 1.215, r ~ 1.1, and a = 5, which produce a mixed phase 
space (see Fig. |2.1| ). The classical results are generated by evolving the initial ensemble 
( |2.25| ) using the mapping ( p.lll) . In Fig. 3^ we compare the time-dependence of the nor- 



malized quantum variance, AL^ = [(L^) — (L)^]//(/ -|- 1), with its classical counterpart, 
AL^ = [(L^)c — (L)^]/|Lp. Circles (diamonds) correspond to the dynamics of an initial 
quantum (classical) state centered at 8{0) = (20°, 40°, 160°, 130°), which is located in the 
connected chaotic zone near one of the unstable fixed points of the classical map. Crosses 
(plus signs) correspond to an initial quantum (classical) state centered on the initial condi- 
tion 9{0) = (5°, 5°, 5°, 5°), which is located in the regular zone near one of the stable fixed 
points. For both initial conditions the quantum and classical results are nearly indistin- 
guishable on the scale of the figure. In the case of the regular initial condition, the quantum 
variance remains narrow over long times and, like its classical counterpart, exhibits a regular 
oscillation. In the case of the chaotic initial condition the quantum variance also exhibits 
a periòdic oscillation but this oscillation is superposed on a very ràpid, approximately ex- 
ponential, growth rate. This exponential growth persists until the variance approaches the 
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Figure 3.1: Growth of normalized quantum and classical variances in a chaotic zone (a) 
and a regular zone (b) of the mixed phase space regime 7 = 1.215 and r ~ 1.1 with / = 154. 
Quantum and classical results are nearly indistinguishable on this scale. In the chaotic case, 
the approximately exponential growth of both variances is governed by a much larger rate, 
= 0.13 (sòlid line), than that predicted from the largest Lyapunov exponent, = 0.04 
(dotted line). 



system size, that is, when AL^ ~ 1 . The initial exponential growth of the quantum variance 
in classically chaotic regimes has been observed previously in several models and appears 
to be a genèric feature of the quantum dynamics; this behaviour of the quantum variance is 
mimicked very accurately by the variance of an initially well-matched classical distribution 



For well-localized states, in the classical case, the exponential growth of the distribution 
variance in chaotic zones is certainly related to the exponential divergence of the underlying 
trajectòries, a property that characterizes classical chaos. To examine this connection we 
compare the observed exponential rate of growth of the widths of the classical (and quantum) 
state with the exponential rate predicted from the classical Lyapunov exponent. For the 
coherent states the initial variance can be calculated exactly, AL^(O) = 1/(1 + 1). Then, 
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assuming exponential growth of this initial variance we get, 

AL^(?7.) ~ y exp(2A^n) for nKtsat-, (3.1) 

where a factor of 2 is included in the exponent since AL^ corresponds to a squared length. 
The dotted line in Fig. 3J^ corresponds to the prediction {"ÒA) with = = 0.04, the 
value of the largest classical Lyapunov exponent. As can be seen from the figure, the actual 
growth rate of the classical (and quantum) variance of the chaotic initial state is significantly 
larger than that predicted using the largest Lyapunov exponent. For comparison purposes 
we also plot a sòlid line in Fig. |3]^ corresponding to (|3.1|) using A^^, = 0.13, which provides 
a much closer approximation to the actual growth rate. We find, for a variety of initial 
conditions in the chaotic zone of this mixed regime, that the actual classical (and quantum) 
variance growth rate is consistently larger than the simple prediction (|3.1| ) using A^ for the 
growth rate. This systematic bias requires some explanation. 



As pointed out in Ref. |46|, the presence of some discrepancy between A^^, and A^ can be 
expected from the fact that the Lyapunov exponent is defined as a geomètric mean of the 
tangent map eigenvalues sampled over the entire connected chaotic zone (corresponding to 
the infinite time limit n ^ oo) whereas the actual growth rate of a given distribution over 
a small number of time-steps will be determined largely by a few eigenvalues of the local 
tangent map. In mixed regimes these local eigenvalues will vary considerably over the phase 
space manifold and the product of a few of these eigenvalues can be quite different from the 
geomètric mean over the entire connected zone. 

However, we find that the actual growth rate is consistently larger than the Lyapunov 
exponent prediction. It is well known that in mixed regimes the remnant KAM tori can 
be 'sticky'; these sticky regions can have a significant decreasing effect on a calculation of 
the Lyapunov exponent. In order to identify an initial condition as chaotic, we specifically 
choose initial states that are concentrated away from these KAM surfaces (regular islands). 
Such initial states will then be exposed mainly to the larger local expansion rates found 
away from these surfaces. This explanation is supported by our observations that, when we 
choose initial conditions closer to these remnant tori, we find that the growth rate of the 
variance is significantly reduced. These variance growth rates are still slightly larger than 
the Lyapunov rate, but this is not surprising since our initial distributions are concentrated 
over a significant fraction of the phase space and the growth of the distribution is probably 
more sensitive to contributions from those trajectòries subject to large eigenvalues away 
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from the KAM boundary than those stuck near the boundary. These explanations are 
further supported by the results of the following section, where we examine a phase space 
regime that is nearly devoid of regular islands. In these regimes we find that the Lyapunov 
exponent serves as a much better approximation to the variance growth rate. 



3.2.2 Regime of Global Chaos 

If we increase the dimensionless couphng strength to 7 = 2.835, with o = 5 and r ~ 1.1 



as before, then the classical flow is predominantly chaotic on the surface V (see Fig. 2.1). 
Under these conditions we expect that genèric initial classical distributions (with non-zero 
support) will spread to cover the full surface V and then quickly relax close to microcanonical 
equilibrium. We find that the initially localised quantum states also exhibit these genèric 
features when the quantum map is governed by parameters that produce these conditions 
classically. 



For the non-autonomous Hamiltonian system (2.11) the total energy is not conserved, 
but the two invariants of motion, and S^, confine the dynamics to the 4-dimensional 
manifold V = S'^xS'^, which is the surface of two spheres. The corresponding microcanonical 



distribution is a constant on this surface, with measure ( 2.12| ), and zero elsewhere. From 



this distribution we can calculate microcanonical equilibrium vàlues for low order moments, 
where, for example, {LJ = {4tt)-^ J^, L^dfi = and {AL^} = {L^} - {L}^ = |L|2. The 
symbols {•} denote a microcanonical average. 

To give a sense of the accuracy of the correspondence between the classical ensemble 



and the quantum dynamics in Fig. 3.2, we show a direct comparison of the dynamics of 
the quantum expectation value (L^) with l = 154 and the classical distribution average 
{Lz)c for an initial coherent state and corresponding classical distribution centered at 9 = 
(45°, 70°, 135°, 70°). To guide the eye in this figure we have drawn lines connecting the 
stroboscopic points of the mapping equations. The quantum expectation value exhibits 
essentially the same dynamics as the classical Liouville average, not only at early times, 
that is, in the initial Ehrenfest regime ||7|, but for times well into the equilibrium regime 
where the classical moment (L^) has relaxed close to the microcanonical equilibrium value 
{Lz} = 0. We have also provided results for a single trajectory launched from the same 
initial condition in order to emphasize the qualitatively distinct behaviour it exhibits. 

In Fig. |3.3| we show the exponential growth of the normalized quantum and classical vari- 



ances on a semilog plot for the same set of parameters and quantum numbers. Numerical 
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Figure 3.2: Comparison of quantum expectation value and corresponding classical average 
{Lz)c in the regime of global chaos 7 = 2.835 and r ~ 1.1 with l = 154 and initial condition 
9o = (45°, 70°, 135°, 70°). The points of the stroboscopic map are connected with lines to 
guide the eye. The quantum expectation value and the Liouville average exhibit esentially 
the same rate of relaxation to microcanonical equilibrium, a behaviour which is qualitatively 
distinct from that of the single trajectory. 



data for (a) correspond to initial condition 9(0) = (20°, 40°, 160°, 130°) and those for (b) cor- 
respond to 0(0) = (45°, 70°, 135°, 70°). As in the mixed regime case, the quantum-classical 
differences are nearly imperceptible on the scale of the figure, and the differences between 
the quantum and classical variance growth rates are many orders of magnitude smaller than 
the small differences in the growth rate arising from the different initial conditions. 



In contrast with the mixed regime case, in this regime of global chaos the prediction (3.1) 
with Xw = = 0.45 now serves as a much better approximation to the exponential growth 
rate of the quantum variance, and associated relaxation rate of the quantum and classical 
states. In this regime the exponent is also much larger than in the mixed regime case 
due to the stronger degree of classical chaos. As a result, the initially localised quantum 
and classical distributions saturate at system size much sooner. 

It is useful to apply ( |3.1j ) to estimate the time-scale at which the quantum (and classical) 
distributions saturate at system size. From the condition AL^ (í^at) — 1 and using (|3.1D we 
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Figure 3.3: Growth of normalized quantum and classical variances in the regime of global 
chaos, 7 = 2.835 and r ~ 1.1 with 1=154, for the two initial conditions cited in the text. 
Quantum-classical differences are nearly imperceptible on this scale. In this regime the 
largest Lyapunov exponent Xl = 0.45 provides a much better estimate of the initial variance 
growth rate. 



obtain, 

tsat-{2XJ-^Hl) (3.2) 

which serves as an estimate of this characteristic time-scale. In the regimes for which the 
full surface V is predominately chaotic, we find that the actual exponential growth rate of 
the width of the quantum state, A^, is well approximated by the largest Lyapunov exponent 
Xl- For a = 5 and r = 1.1, the approximation ~ A^ holds for coupling strengths 7 > 2, 
for which more than 99% of the surface V is covered by one connected chaotic zone (see 



Fig-iD 



By comparing the quantum probability distribution to its classical counterpart, we can 
learn much more about the relaxation properties of the quantum dynamics. In order to 
compare each mi value of the quantum distribution, Pz{mi), with a corresponding piece of 
the continuous classical marginal probability distribution, 

PciL,) = í íí dS,d(l),d(l)i pci9s, <Ps, Oi, (Pi), (3.3) 
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Figure 3.4: Initial probability distributions for for 6^(0) = (45°, 70°, 135°, 70°) with 
/ = 154. The quantum and classical distributions are indistinguishable on the scale of the 
figure. 



we discretize the latter into 2j + 1 bins of width h = 1. This procedure produces a discrete 
classical probabiüty distribution P^{mi) that prescribes the probabihty of finding the spin 
component in the interval [m/ + 1/2, mi — 1/2] along the z-axis. 

To illustrate the time-development of these distributions we compare the quantum and 
classical probability distributions for three successive vàlues of the kick number n, using the 



same quantum numbers and initial condition as in Fig. |3.2| . In Fig. 3.4 the initial quantum 
and classical states are both well-localised and nearly indistinguishable on the scale of the 
figure. At time n = 6 ~ tgat, shown in Fig. both distributions have grown to fill 
the accessible phase space. It is at this time that the most significant quantum-classical 
discrepàncies appear. 

For times greater than tgat, however, these emergent quantum-classical discrepencies do 
not continue to grow, since both distributions begin relaxing towards equilibrium distribu- 
tions. Since the dynamics are confined to a compact phase space, and in this parameter 
regime the remnant KAM tori fill a negligibly small fraction of the kinematicaly accessi- 
ble phase space, we might expect the classical equilibrium distribution to be very close to 
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6 in the regime of global 

2.835 and r ~ 1.1. Both the quantum and classical distribution have spread to 
system dimension and exhibit their largest differences on this saturation time-scale. 



Figure 3.5: Same as Fig. 3.4 but the states have evolved to n 
chaos 7 



the microcanonical distribution. Indeed such relaxation close to microcanonical equilibrium 
is apparent for both the quantum and the classical distribution at very early times, as 
demonstrated in Fig. corresponding to n = 15. 

Thus the signature of a classically hyperbolic flow, namely, the exponential relaxation of 
an arbitrary distribution (with non-zero measure) to microcanonical equilibrium [^0[, holds 
to good approximation in this model in a regime of global chaos. More suprisingly, this 
classical signature is manifest also in the dynamics of the quantum distribution. In the 



quantum case, however, as can be seen in Fig. 3.6, the probability distribution is subject to 
small irreducible time-dependent fluctuations about the classical equilibrium. We examine 
these quantum fluctuations in detail Chapter 4. 
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Figure 3.6: Same as Fig. 3.5, but for n = 15. Both quantum and classical distributions 



have relaxed close to the microcanonical equilibrium. 



3.3 Time-Domain Characteristics of Quantum-Classical 
Differences 

We consider the time dependence of quantum-classical differences defined along the z-axis 
of the spin L, 

6L,{n) = \{L,{n))-{L,{n))cl (3.4) 



at the stroboscopic times t = n. In Fig. |3.7| we compare the time-dependence of 5Lz{n) 
on a semi-log plot for a chaotic state (filled circles), with 0(0) = (20°, 40°, 160°, 130°), and 
a regular state (open circles), 0{tí) = (5°, 5°, 5°, 5°), evolved using the same mixed regime 



parameters (7 = 1.215 and r ~ 1.1) and quantum numbers (/ = 154) as in Fig. 3.1 



We are interested in the behaviour of the upper envelope of the data in Fig. |3.7| . For the 
regular case, the upper envelope of the quantum-classical differences grows very slowly, as 
some polynomial function of time. For the chaotic case, on the other hand, at early times 
the difference measure (|3.4| ) grows exponentially until saturation around n = 15, which is 
well before reaching system dimension, |L| ~ / = 154. After this time, which we denote 
í*, the quantum-classical differences exhibit no definite growth, and fluctuate about the 
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Figure 3.7: Time-dependence of quantum-classical differences in a regular zone (open 
circles) and a chaotic zone (filled circles) of mixed regime (7 = 1.215 and r ~ 1.1) with 
/ = 154. For the chaotic state, 6Lz = \{Lz) — {Lz)c\ is contrasted with the Ehrenfest 
difference \ {Lz) — Lz\ between the quantum expectation value and a single trajectory (plus 
signs), which grows until saturation at system dimension. The sòlid line corresponds to 
( |3.5D using Xqc = 0.43. The horizontal lines indicate two different vàlues of the difference 
tolerance p which may be used to determine the break-time; for p = 0.1 (dotted line) tf, 
occurs on a logarithmic time-scale, but for p = 15.4 (sparse dotted line) is not defined 
over numerically accessible time-scales. 



equilibrium value óL^ ~ 1 ^ |L|. In Fig. 3.7 we also include data for the time-dependence 
of the Ehrenfest difference | (L^) — , which is defined as the difference between the quantum 
expectation value and the dynamical variable of a single trajectory initially centered on the 
quantum state. In contrast to ÓL^, the ràpid growth of the Ehrenfest difference continues 
until saturation at the system dimension. 



In Fig. |3.8| we compare the time-dependence of the quantum-classical differences in 
the case of the chaotic initial condition 6{0) = (20°, 40°, 160°, 130°) for quantum numbers 



/ = 22 (filled circles) and / = 220 (open circles), using the same parameters as in Fig. 3.7 . 
This demonstrates the remarkable fact that the exponential growth terminates when the 
difference measure reaches an essentially fixed magnitude (óL^ ~ 1 as for the case / = 154), 
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Figure 3.8: Growth of the quantum-classical difference óL^ in the chaotic zone of a mixed 
regime, 7 = 1.215 and r ~ 1.1, with / = 22 (filled circles) and l = 220 (open circles). For 
l = 220 the exponential growth rate (^) is plotted using the classical Lyapunov exponent, 
Xl = 0.04 (sparse dotted hne), and for both / vàlues ( ^.51) is plotted using the exponent 
Agc = 0.43 (sòlid line for l = 22, dotted hne for / = 220), which is obtained from a fit of 
( |3.6D to the corresponding break-time data in Fig. ^.10| . 



although the system dimension differs by an order of magnitude in the two cases. 



In Fig. 3.£ we consider the growth of the quantum-classical difference measure 5Lz{n) in 
a regime of global chaos, for / = 154, and using the same set of parameters as those examined 
in Fig. |3.3| (7 = 2.835 and r ~ 1.1). Again the upper envelope of the difference measure 
5Lz[n) exhibits exponential growth at early times, though in this regime of global chaos 
the exponential growth persists only for a very short duration before saturation at t* ~ 6. 
The initial condition 0(0) = (20°, 40°, 160°, 130°) is a typical case (fiUed circles), where, as 
seen for the mixed regime parameters, the magnitude of the difference at the end of the 
exponential growth phase saturates at the value 5Lz{t*) ~ 1, which does not scale with the 
system dimension (see Fig. 3.11). The initial condition 9(0) = (45°, 70°, 135°, 70°) (open 
circles) leads to an anomolously large deviation at the end of the exponential growth phase, 
6Lz{t*) ~ 10, though still small relative to the system dimension |L| ~ 154. This deviation 
is transient however, and at later times the magnitude of quantum-classical differences 
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Figure 3.9: Growth of quantum-classical differences in the regime of global chaos, 7 = 2.835 
and r ~ 1.1, with / = 154, for the two initial conditions cited in text. The exponential growth 
rate ( |3.5[ ) is plotted using the classical Lyapunov exponent, Xl = 0.45 (dotted line), and the 
exponent Xqc = 1.1 (sòlid line), which is obtained from a fit of (|3.6| ) to the corresponding 
break-time data in Fig. |3.10| . 



fluctuates about the equilibrium value 6Lz ~ 1. The quantum-classical differences are a 
factor of l/l smaller than typical differences between the quantum expectation value and 



the single trajectory, which are of order system dimension (see Fig. 3.2) as in the mixed 
regime case. 

In all cases where the initial quantum and classical states are launched from a chaotic 
zone we find that the initial time-dependence of quantum-classical differences compares 
favorably with the exponential growth ansatz. 



tó^(n) ~ ■^exp(Agcn) for n<t*, (3.5) 
8í 



where the exponent Xqc is a new exponent subject to numerical measurement |11]. Since 
contributions from the initial differences in other mismatched moments will generally mix 
under the dynamical flow, it is appropriate to consider an effective initial difference for the 
prefactor in ( |3.5| ). The prefactor 1/81 is obtained by accounting for the initial contributions 
from the 3 cartesian components, [^^^^(O) + S^Ly{0) + Ó^L^iO)^/'^ = 1/81. 
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We are interested in whether the Lyapunov exponent Al is a good approximation to Xgc- 
In Fig. U we plot with A^c = = 0.04 (dotted line) for / = 220. Clearly the largest 
Lyapunov exponent severly underestimates the exponential growth rate of the quantum- 
classical differences, in this case by more than an order of magnitude. The growth rate of 
the state width, A^ = 0.13 , is also several times smaller than the initial growth rate of the 
quantum-classical differences. In the case of Fig. |3.9| , corresponding to a regime of global 
chaos with a much larger Lyapunov exponent, we plot ( |3.5| ) with Agc = Al = 0.45 (dotted 
line), demonstrating that, in this regime too the largest Lyapunov exponent underestimates 
the initial growth rate of the quantum-classical difference measure 6Lz{n). 

We also find, from inspection of our results, that the time t* at which the exponential 
growth ( ^.51) terminates can be estimated from tgat, the time-scale on which the distributions 



saturate at or near system size (3.2). In the case of the chaotic initial condition of Fig. 3.1 



for which 7 = 1.215, visual inspection of the figure suggests that tgat — 18. This should be 
compared with Fig. \i.7[ where the exponential growth of 5Lz{n) ends rather abruptly at 
t* ~ 15. In Fig. |3.3| , corresponding to a regime of global chaos (7 = 2.835), the variance 
growth saturates much earlier, around tgat — 6 for both initial conditions. From Fig. p.9| it is 
apparent that in this regime t* ~ 6. As we increase 7 further, we find that the exponential 
growth phase of quantum-classical differences 6Lz{n) is shortened, lasting only until the 
corresponding quantum and classical distributions saturate at system size. For 7 ~ 12, with 
Al — 1.65, the chaos is sufficiently strong that the initial coherent state for l = 154 spreads 
to cover V within a single time-step. Similarly the initial difference measure 6Lz{0) ~ 0.001 
grows to the magnitude 6Lz{l) — 1 within a single time-step and subsequently fluctuates 
about that equilibrium value. We have also inspected the variation of t* with the quantum 



numbers and found it to be consistent with the logarithmic dependence of tsat in (3.2). 



3.4 Correspondence Scaling in the Macroscopic Limit 

We have assumed in ( |3.5D that the exponent Agc is independent of the quantum numbers. 
A convenient way of confirming this, and also estimating the numerical value of Xgc, is by 
means of a break-time measure. The break-time is determined from the time tij{l,p) at 
which quantum-classical differences exceed some tolerance p, which does not scale with the 
system size for each quantum model. The classical parameters and initial condition are held 
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fixed. Setting óLz{ti)) = p in (^), we obtain tt, in terms of p, l and Agc, 

h ~ ln(8 p l) provided p < 0(1). (3.6) 

The restriction p < 0(1), which plays a crucial role in limiting the robustness of the break- 
time measure (|3.6| ), is explained and motivated further below. 

The exphcit form we have obtained for the argument of the logarithm in (3.6) is a 
direct result of our estimate that the initial quantum-classical differences arising from the 
Cartesian components of the spin provide the dominant contribution to the prefactor of the 
exponential growth ansatz (|3.5| ). Differences in the mismatched higher order moments, as 
weh as intrinsic differences between the quantum dynamics and classical dynamics, may also 
contribute to this effective prefactor. We have checked that the initial value 5Lz{tí) ~ 1/8/ is 
an adequate estimate by comparing the intercept of the quantum-classical data on a semilog 



plot with the prefactor of (|3.5D for a variety of / vàlues (see e.g. Fig. 3.^ 



In Fig. we examine the scaling of the break-time for / vàlues ranging from 11 to 
220 and with fixed tolerance p = 0.1. The break-time can assume only the integer vàlues 
t = n and thus the data exhibits a step-wise behaviour. For the mixed regime parameters, 
7 = 1.215 and r ~ 1.1 (filled circles), with initial condition ^(0) = (20", 40°, 160°, 130°), a 
non-linear least squares fit to ( |3.6| ) gives Xqc = 0.43. This fit result is plotted in the figure as 
a sòlid line. The close agreement between the data and the fit provides good evidence that 
the quantum-classical exponent Xqc is independent of the quantum numbers. To check this 



result against the time-dependent 5Lz{n) data, we have plotted the exponential curve (3.5) 
with Açc = 0.43 in Fig. |3.7| using a sòlid line and in Fig. using a sòlid line for / = 22 and 
a dotted line for / = 220. The exponent obtained from fitting ( |3.6D serves as an excel·lent 
approximation to the initial exponential growth ( p.5D of the quantum-classical differences 
in each case. 

In Fig. |3.10| we also plot break-time results for the global chaos case 7 = 2.835 and 
r ~ 1.1 (open circles) with initial condition 9{tí) = (45°, 70°, 135°, 70°). In this regime the 
quantum-classical differences grow much more rapidly and, consequently, the break-time 
is very short and remains nearly constant over this range of computationally accessible 
quantum numbers. Due to this limited variation, in this regime we can not confirm ( |3.6D , 
although the data are consistent with the predicted logarithmic dependence on /. Moreover, 
the break-time results provide an effective method for estimating Xqc if we assume that 
(|3.6[) holds. The same fit procedure as detailed above yields the quantum-classical exponent 
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Figure 3.10: Scaling of the break-time using tolerance p = 0.1 as a function of increasing 
quantum number for the mixed regime parameters 7 = 1.215 and r ~ 1.1 with 9{0) = 
(20°, 40°, 160°, 130°) (filled circles) and for the global chaos parameters 7 = 2.835 and 
r ~ 1.1 with (9(0) = (45°, 70°, 135°, 70°) (open circles). We also plot the results of fits 
to the log rule (|3.6|), which produced exponents Xqc = 0.43 for 7 = 1.215 and Xgc = 1.1 for 
7 = 2.835. 



Xqc = 1.1. This fit result is plotted in Fig. |3.10| as a sòlid line. More importantly, the 
exponential curve ( |3.5D , plotted with fit result Xqc = 1.1, can be seen to provide very good 



agreement with the initial growth rate of Fig. 3.9 for either initial condition, as expected. 

In the mixed regime (7 = 1.215), the quantum-classical exponent Xqc = 0.43 is an order 
of magnitude greater than the largest Lyapunov exponent Al = 0.04 and about three times 
larger than the growth rate of the width Xw = 0.13. In the regime of global chaos (7 = 2.835) 
the quantum-classical exponent Xqc = 1.1 is a little more than twice as large as the largest 
Lyapunov exponent = 0.45. 

The condition p < 0(1) is a very restrictive limitation on the domain of application 
of the log break-time ( |3.6| ) and it is worthwhile to explain its significance. In the mixed 
regime case of Fig. |3.7| , with / = 154, we have plotted the tolerance vàlues p = 0.1 (dotted 
line) and p = 15.4 (sparse dotted line). The tolerance p = 0.1 is exceeded at í = 11, while 
the quantum-classical differences are still growing exponentially, leading to a log break-time 
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for this tolerance value. For the tolerance p = 15.4 <C |L|, on the other hand, the break- 
time does not occur on a measurable time-scale, whereas according to the logarithmic rule 



(3.6), with l = 154 and Xqc = 0.43, we should expect a rather short break-time ~ 23. 
Consequently the break-time ( |3.6D , applied to delimiting the end of the Liouvihe regime, is 
not a robust measure of quantum-classical correspondence. 

Our definition of the break-time (^]^) requires holding the tolerance p fixed as the system 
size increases (and not as a fraction of the system dimension as in [^]) when comparing 
systems with different quantum numbers. Had we chosen to compare systems using a 
fixed relative tolerance, /, then the break-time would be of the form íf, ~ A^^^ ln(8 / 
and subject to the restriction / < 0{l/l). Since / — > in the classical limit, this form 
emphasizes that the log break-time applies only to differences that are vanishing fraction of 
the system dimension in that limit. 

Although we have provided numerical evidence (in Fig. |3.8| ) of one mixed regime case in 
which the largest quantum-classical differences occuring at the end of the exponential growth 
period remain essentially constant for varying quantum numbers, 6Lz{t*) ~ C(l), we find 
that this behaviour represents the typical case for all parameters and initial conditions 
which produce chaos classically. To demonstrate this behaviour we consider the scaling 
(with increasing quantum numbers) of the maximum vàlues attained by 5Lz{n) over the 
first 200 kicks, (5L™"^. Since t* ^ 200 over the range of l vàlues examined, the quantity 
^jmax -g ^ rigorous upper bound for 6Lz{t*). 



In Fig. 3.11 we compare JL™"^ for the two initial conditions of Fig. |3.9| and using the 
global chaos parameters (7 = 2.835, r ~ 1.1). The filled circles in Fig. 3.11| correspond to 
the initial condition 9(0) = (20°, 40°, 160°, 130°). As in the mixed regime, the maximum 
deviations exhibit little or no scaling with increasing quantum number. This is the typical 
behaviour that we have observed for a variety of different initial conditions and parameter 
vàlues. These results motivate the genèric rule. 



' max 



SLAt') <«L^^ 0(1/0. (3.7) 



Vi(ÏTT) - ■M+ï) 

Thus the magnitude of quantum-classical differences reached at the end of the exponential 
growth regime, expressed as a fraction of the system dimension, approaches zero in the 
classical limit. 

However, for a few combinations of parameters and initial conditions we do observe a 
'transient' discrepancy peak occuring at í ~ í* that exceeds 0{1). This peak is quickly 
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Figure 3.11: Maximum quantum-classical difference occuring over the first 200 kicks in the 
regime of global chaos (7 = 2.835, r ~ 1.1) plotted against increasing quantum number. 
These maximum vàlues provide an upper bound on 6Lz{t*) for each /. The data correspond- 
ing to the initial condition 9{0) = (20°, 40°, 160°, 130°) (filled circles) represent a typical case 
in which the maximum quantum-classical differences do not vary significantly with l. The 
large deviations observed for the initial condition 9(0) = (45°, 70°, 135°, 70°) (open circles) 
are an exceptional case, with maximum differences growing rapidly for small quantum num- 
bers but tending asymptotically toward independence of /. These curves provide an upper 
bound on the tolerance vàlues p for which the break-time measure scales logarithmically 
with /. 



smoothed away by the subsequent relaxation of the quantum and classical distributions. 
This peak is apparent in Fig. |3.9| (open circles), corresponding to the most conspicuous 
case that we have identified. This case is apparent as a small deviation in the normalized 
data of Fig. ^.2| . The scaling of the magnitude of this peak with increasing / is plotted 
with open circles in Fig. |3.11 . The magnitude of the peak initially increases rapidly but 
appears to become asymptotically independent of /. The other case that we have observed 
occurs for the classical parameters 7 = 2.025, with r ~ 1.1 and a = 5, and with initial 
condition 9(0) = (20°, 40°, 160°, 130°). We do not understand the mechanism leading to 
such transient peaks, although they are of considerable interest since they provide the most 
prominent examples of quantum-classical discrepancy that we have observed. 
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3.5 Discussion 

In this study of a non-integrable model of two interacting spins we have characterized the 
correspondence between quantum expectation vàlues and classical ensemble averages for 
intially localised states. We have demonstrated that in chaotic states the quantum-classical 
differences initially grow exponentially with an exponent Xqc that is consistently larger than 
the largest Lyapunov exponent. In a study of the moments of the Henon-Heiles system, 
Ballentine and McRae ||ïï|, have also shown that quantum-classical differences in chaotic 
states grow at an exponential rate with an exponent larger than the largest Lyapunov 
exponent. This exponential behaviour appears to be a genèric feature of the short-time 
dynamics of quantum-classical differences in chaotic states. 

Since we have studied a spin system, we have been able to solve the quantum problem 
without truncation of the Hilbert space, subject only to numerical roundoff, and thus we are 
able to observe the dynamics of the quantum-classical differences well beyond the Ehrenfest 
regime. We have shown that the exponential growth phase of the quantum-classical differ- 
ences terminates well before these differences have reached system dimension. We find that 
the time-scale at which this occurs can be estimated from the time-scale at which the dis- 
tribution widths approach the system dimension, tsat — (2AtL,)^^ ln(/) for initial minimum 
uncertainty states. Due to the close correspondence in the growth rates of the quantum 
and classical distributions, this time-scale can be estimated from the classical physics alone. 
This is useful because the computational complexity of the problem does not grow with the 
system action in the classical case. Moreover, we find. tliB-t tliG exponent A^^ can be ap- 
proximated by the largest Lyapunov exponent when the kinematic surface is predominantly 
chaotic. 

We have demonstrated that the exponent Xqc governing the initial growth rate of the 
quantum-classical differences is independent of the quantum numbers, and that the effective 
prefactor to this exponential growth decreases as l/l. These results imply that a log break- 
time rule (^^) delimits the dynamical regime of Liouville correspondence. However, the 
exponential growth of quantum-classical differences persists only for short times and small 
differences, and thus this log break-time rule applies only in a similarly restricted domain. 
In particular, we have found that the magnitude of the differences occuring at the end of 
the initial exponential growth phase does not scale with the system dimension. A typical 
magnitude for these differences, relative to the system dimension, is 0{l/l). Therefore, log(/) 
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break-time rules characterizing the end of the Liouville regime are not robust, since they 
apply to quantum-classical differences only in a restricted domain, i.e. to relative differences 
that are smaller than 0{l/l). 

This restricted domain effect does not arise for the better known log break-time rules 



describing the end of the Ehrenfest regime 15, 54 1. The Ehrenfest log break-time remains 
robust for arbitrarily large tolerances since the corresponding differences grow roughly expo- 
nentially until saturation at the system dimension p^ , pTt - Consequently, a log(/) break-time 
indeed implies a breakdown of Ehrenfest correspondence. However, the logarithmic break- 
time rule characterizing the end of the Liouville regime does not imply a breakdown of 
Liouville correspondence because it does not apply to the observation of quantum-classical 
discrepàncies larger than 0(1//). The appearance of residual 0(1//) quantum-classical dis- 
crepàncies in the description of a macroscopic body is, of course, consistent with quantum 
mechanics having a proper classical limit. 

We have found, however, that for certain exceptional combinations of parameters and 
initial conditions there are relative quantum-classical differences occuring at the end of the 
exponential growth phase that can be larger than 0(1//), though still much smaller than 
the system dimension. In absolute terms, these transient peaks seem to grow with the 
system dimension for small quantum numbers but become asymptotically independent of 
the system dimension for larger quantum numbers. Therefore, even in these least favorable 
cases, the fractional differences between quantum and classical dynamics approach zero in 
the limit / — > oo. This vanishing of fractional differences is sufficient to ensure a classical 
limit for our model. 

Finally, contrary to the results found in the present model, it has been suggested that 
a log break-time delimiting the Liouville regime implies that certain isolated macroscopic 
bodies in chaotic motion should exhibit non-classical behaviour on observable time scales. 
However, since such non-classical behaviour is not observed in the chaotic motion of macro- 
scopic bodies, it is argued that the observed classical behaviour emerges from quantum 
mechanics only when the quantum description is expanded to include interactions with the 



many degrees-of-freedom of the ubiquitous environment ||97| , 99 1. (This effect, called de- 
coherence, rapidly evolves a pure system state into a mixture that is essentially devoid of 
non-classical properties.) However, in our model the classical behaviour emerges in the 
macroscopic limit of an isolated few degree-of-freedom quantum system that is described by 
a pure state and subject only to unitary evolution. Quantum-classical correspondence at 
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both early and late times arises in spite of the log break-time because this break-time rule 
applies only when the quantum-classical difference threshold is chosen smaller than 0{h). 
In this sense we find that the decoherence effects of the environment are not necessary for 
correspondence in the macroscopic hmit. Of course the effect of decoherence may be exper- 
imentally significant in the quantum and mesoscopic domains, but it is not required as a 
matter of principle to ensure a classical hmit. 



Chapter 4 

Correspondence for the Probability 
Distributions 



This chapter is taken from Emerson and Ballentine |36| 



4.1 Introduction 

The study of chaos in quantum dynamics has led to differing views on the conditions required 
for demonstrating quantum-classical correspondence [Q, |5^. Moreover, the criteria by which 



this correspondence should be measured have also been a subject of some controversy 23, 
p7| . While much of the earher work on this tòpic is concerned with characterizing the degree 
of correspondence between quantum expectation vàlues and classical dynamical variables 
[ p!5| , p8| , the more recent approach is to focus on differences between the properties of 
quantum states and associated classical phase space densities evolved according to Liouville's 
equation 0, |5|, ||, ^, |ïï|, ||. 



Several authors have examined quantum-classical correspondence by considering the ef- 



fects of interactions with a stochastic environment |5^, 52], a process sometimes called 
decoherence. While this process may improve the degree of quantum-classical correspon- 
dence for fixed quantum numbers, it has been further suggested that the limit of large 
quantum numbers is inadequate for correspondence, and that decoherence must be taken 
into account to generate classical appearances from quantum theory; this view has been 
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argued to apply even in the case of macroscopic bodies that are described initially by well- 



localised states, provided their classical motion is chaotic [97, In this chapter we 

examine how the degree of correspondence with Liouville dynamics scales specifically in the 
hmit of large quantum numbers. This "classical limit" is distinct from a "thermodynamic 
limit", that is, a limit involving many quantum numbers. 

The degree of Liouville correspondence has been characterized previously by studying the 
differences between the means and variances of the dynamical variables ||7, 57, |7^, 11, |3^, 12]. 
This involves a comparison of quantum expectation vàlues and classical ensemble averages. 
However, these low-order moments give only crude Information about the differences be- 
tween the quantum and classical states. Specifically, the quantum state may exhibit coarse 
structure which differs significantly from the classical state although the means and vari- 
ances (for some simple observabes) are nearly the same for the quantum and classical states. 
Moreover, much of the previous work was concerned with correspondence at early times, 
or more precisely, in the Ehrenfest regime when the states are narrow compared to system 
dimensions ^ . 

Another approach is to identify quantum-classical differences with differences between 



the Wigner quasi-distribution and the classical phase space density |56|. This approach 
is objectionable because the Wigner quasi-distribution may take on negative vàlues and 
therefore may not be interpreted as a "classically observable" phase space distribution. It is 
possible to consider instead smoothed quantum phase space distributions, but in this case 
the residual quantum-classical differences still do not have clear experimental significance. 

In this chapter, we characterize the degree of quantum-classical correspondence by com- 
paring quantum probability distributions for dynamical variables with the corresponding 
classical marginal distributions for these dynamical variables. These are well-defined classi- 
cal observables that describe the distribution of outcomes upon measurement of the given 
dynamical variable. We are interested in the differences that arise on a fine scale and 
therefore characterize the typical quantum-classical deviations that arise in bins of width h. 

The dynamics are generated by the model of interacting spins described in Chapter 2. 
The Hilbert space is finite dimensional so no artificial truncation of the state is required. 
The quantum time-evolution is unitary and the classical motion is volume-preserving (sym- 
plectic). In the case of classically chaotic motion, we follow initially localised states until 
they have evolved well beyond the relaxation time-scale of the classical density. Throughout 
the chapter we emphasize that the quantum signatures of chaos that appear in the quantum 
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distributions are the same as those that appear in the marginal classical distributions. In 
particular, the quantum relaxation rates can be accurately estimated from the Liouville 
dynamics of an approximately matching initial phase space density. This purely classi- 
cal approximation is surprisingly accurate even for small quantum numbers, but may be 
most useful for the theoretical description of mesoscopic systems since the purely classical 
calculations do not scale with the quantum numbers. 

The quantum and classical probability distributions remain close even after the states 
have spread to the system dimension. Specifically, in mixed regimes, the quantum distribu- 
tions exhibit an equilibrium shape that reflects the detalls of the classical KAM surfaces. 
When the classical manifold is predominantly chaotic, the quantum and classical states relax 
close to the microcanonical state. However, in both of these chaotic regimes the equilib- 
rium quantum distributions exhibit characteristic fluctuations away from the classical ones. 
We demonstrate that the Standard deviation of these quantum-classical differences becomes 
vanishingly small in the classical limit, J /h oc, where JT" is a characteristic system action. 



4.2 Dynamical Behaviour of Probability Distributions 

In the case of a classical mixing system, initial densities with non-zero measure are expected 
to spread in an increasingly uniform manner throughout the accessible phase space. The 
term uniform is meant to apply specifically in a coarse-grained sense. For some simple maps, 
such as the baker's map, it is possible to show that this rate of relaxation to the equilibrium 
configuration occurs exponentially with time pO|]. 



The spin map we consider ( 2.11 ) is not mixing on the accessible classical manifold V, 
but has mixed dynamics: depending on the system parameters, the surface V can generally 
be decomposed into regions of regular dynamics and a connected region of chaotic dynamics 
[p5| . In parameter regimes that are predominantly chaotic, we expect behaviour on V that 
approximates that of a mixing system. In particular, initially localised Liouville densities 
should relax close towards the microcanonical measure at an exponential rate, on average. 
In this section we demonstrate that these signatures of chaos are exhibited also by the 
quantum dynamics. Most striking is the degree of similarity between the quantum and 
classical behaviours even in regimes with classically mixed dynamics. 

We are interested in the behaviour of quantum probability distributions that are as- 
sociated with measurements of classical dynamical variables. The quantum probability 
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distribution associated with the classical observable is given by, 

PlMi) = {Hn)\ Ri,rrn Wn)) = Tr [ \l,mi){l,mi\p'^'\n) ] , (4.1) 

where, 

Ri,mi = l.s'^\l,mi){l,mi\ (4.2) 
is a projection operator onto the eigenstates of L^, and 

/)(n) = TVW[|V(n))(V^(n)|], (4.3) 

is the reduced state operator for the spin L at time n and Tr^*-* denotes a trace over the 
factor space Tís- We have written out the exphcit expression (0) to emphasize that the 
probabihty of obtaining each mi value is associated with a projector onto a subspace of the 
factor space Tíi. 

For reasons related to this fact (which we will make clear in later sections), we are 
also interested in examining the probabihty distributions associated with components of the 
total angular momentum J = S + L. The probability of obtaining a given rrij value upon 
measurement of Jz is given by, 

PjÀ^j) = XlKV'('^)l·s,í,"^s>"^i - (4.4) 



where |s, Z, m^, m-,- — rris) is an element of the orthonormal basis (p.2|). The probability 
Pj^{mj) is associated with a projector onto a subspace of the full Hilbert space Tí. The 
dimension of each subspace is given by the number of pairs {nisjmi) that yield a given value 
of rrij = rUs + mi . 

The classical probability distributions associated with dynamical variables are obtained 
by partial integration over the accessible phase space. In the case of L^, the continuous 
marginal distribution is given by, 

P{Lz) = í í í dSzdcl)sdct)i p,{Sz,cl)s,Lz,cl)i), (4.5) 



where for notational convenience we have suppressed reference to the time-dependence. The 
marginal probability distribution for the total spin component Jz is obtained by integration 
subject to the constraint Sz + Lz = Jz, 

PiJz) =1111 dSzd<l)sdLzd(t>i pc{Sz,^s,Lz,^i) 6{Sz + Lz - Jz). (4.6) 
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These classical distributions are continuous, though their quantum counter-parts are in- 
trinsically discrete. To construct a meaningful quantum-classical comparison it is useful to 
discretize the classical distrbutions by integrating the continuous probabilities over intervals 
of width h = 1 centered on the quantum eigenvalues. In the case of the component L^, the 
quantum probability PL^{mi) is then associated with the classical probability of finding Lz 
in the interval [rrii — 1/2, m/ + 1/2]. This is given by 

fmi+l/2 

PLimi)= P{Lz). (4.7) 

Similarly, in the case of J^, we compare each quantum Pj^{raj) with the discrete classical 
probability, 

Pl{mj)= P{Jz). (4.8) 

In the following discussion of the numerical results we will emphasize that, for chaotic 
states, the steady-state shape of the quantum and classical distributions should be com- 
pared with the marginal distributions derived from the microcanonical state. Our model is 
non-autonomous, but the spin magnitudes are conserved. The appropriate classical micro- 
canonical measure is a constant on the accessible manifold V = S"^ x S'^. This follows from 
the usual equilibrium hypothesis that all accessible microstates are equiprobable, where 
equiprobability is defined with respect to the invariant measure ( 2.12| ). This microcanonical 



density projected onto the L^-axis produces the discrete, flat distribution, 

Prf{mi) = {2l + l)-\ (4.9) 

However, projected along J^, the microcanonical distribution is not flat, but has a tent- 
shape, 

Primi) = l + '+^-J'^^l for \mA>l-s 
■^^ ^ (2s + 1)(2/ + 1) ' - 

= 2l\-l \mj\<l-s. (4.10) 

In quantum mechanics, the equiprobability hypothesis implies that the appropriate micro- 
canonical state is an equal-weight mixture. This microcanonical state, sometimes called a 
random state, is proportional to the identity in the full Hilbert space 7i = Tíg ® 'Hi- It 
produces the same projected microcanonical distributions, i.e. ( |4.9| ) for Lz and ( 4.10| ) for 
J2, as the classical microcanonical state. 
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Figure 4.1: Quantum and classical probability distributions for with l = 154 in chaotic 
zone of mixed regime (7 = 1.215, r = 1.1, a = 5). The dots are visible because they are 
shifted to the right by half of their width. The figure on the left is the initial state ( n = 0) 
and that on the right is at time-step n = 6. 



4.2.1 Mixed Regime Chaos 

We consider first a classical parameter regime (7 = 1.215, r = 1.1, and a = 5) for which the 
kinematically accessible phase space V is highly mixed. The chaotic region appears to be 
connected (all chaotic initial conditions have the same largest Lyapunov exponent A^, = 0.04) 
and covers about half of the kinematic surface. A projection of only the chaotic initial 
conditions onto the plane spanned by Sz and reveals large regular islands surrounding 
the stable paral·lel fixed points (itS^, itL^), with chaotic regions spreading out from the 
unstable anti-paral·lel fixed points {zizSz,^Lz)- A similar projection of the regular initial 
conditions shows points not only clustered about the paral·lel fixed points but also spread 
along the axis Sz = Lz- 

We now consider the time-evolution of quantum and classical states concentrated in the 
chaotic zone near one of the unstable anti-parallel fixed points, with initial centroids directed 
along 0(0) = (20°, 40°, 160°, 130°). The quantum dynamics are calculated using quantum 
numbers s = 140 and / = 154. As shown in Fig. 4.1, at early times both the quantum 



distribution Pi^[mi) (sòlid line) and the corresponding classical distribution Pf^^[mi) (dots) 
remain well-localised. Their initial differences are not distinguishable on the scale of the 
figures. (The dots are shifted to the right by half of their width.) By time-step n = 20 both 
quantum and classical distributions have broadened to the system dimension and begin to 
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Figure 4.2: Same as Fig. 4.1 but for time-steps n = 99 on the left and n = 100 on the 
right. Both quantum and classical distributions have reached the system dimension and are 
relaxing towards equilibrium. 



exhibit noticeable differences. As shown in Fig. around n = 100 the distributions have 
begun to settle close to an equihbrium shape. In Fig. 4.3 the successive time steps n = 199 



and n = 200 show that, although both the quantum and classical distributions have relaxed 
very close to the same equilibrium distribution, the quantum distribution exhibits rapidly 
oscillating fluctuations about the classical steady-state. 

Both the quantum and classical equilibrium distributions (projected along Lz) show 
significant deviation from the microcanonical distribution ( |4.9| ). This is also true of the 
distribution projected along L^, which has a diíFerent non-uniform equilibrium distribu- 



tion than that observed when projecting onto Lz (see the left box of Fig. 4.4). Uniform 
marginal distributions would be expected if the classical mapping was mixing, in which case 
arbitrary initial densities (with non-vanishing measure) would relax to the microcanonical 
distribution. Since the accessible kinematic surface has large KAM surfaces in this parame- 
ter regime, the coarse-grained classical equilibrium distributions are not expected to be flat. 
An unexpected feature of the results is the observation that the shape of the equilibrium 
quantum distributions so accurately reflects the detalls of the KAM structure in the classical 
phase space. This feature is most striking in the case of the distributions projected along 
Jz (see the right box in Fig. iA). The steady-state quantum and classical probability dis- 
tributions Pj^(mj) and Pj^{mj) are both sharply peaked about rrij = 0. This equilibrium 
shape is much more sharply peaked than the tent-shape of the projected microcanonical 



distribution, Pj^'^(mj), given by ( 4.1C| ) and also plotted in the right box in Fig. iÀ. The 
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Figure 4.3: Same as Fig. 4.1 but for n = 199 on the left and n = 200 on the right. The 
quantum distribution is fluctuating about a classical steady-state. 



important point is that the additional localization of the quantum distribution can be under- 
stood from a Standard fixed-point analysis of the classical map [^ ]: the presence of KAM 
surfaces arising due to the stabihty of the parallel fixed points prevents the chaotic classical 
spins from aligning in parallel along the z-axis. Most remar kably, we find that the steady- 
state quantum distributions accurately reproduce this parameter-dependent structure of the 
mixed classical phase space even for much smaller quantum numbers. We examine how the 



accuracy of this correspondence scales with the quantum numbers in section 4.5 



4.2.2 Regime of Global Chaos 

If we hold o = 5 and r = 1.1 fixed and increase the coupling strength to the value 7 = 
2.835, then all four of the fixed points mentioned above become unstable []35| . Under these 
conditions less than 0.1% of the surface V is covered with regular islands; the remainder 
of the surface produces a connected chaotic zone with largest Lyapunov exponent = 
0.45. We will sometimes refer to this parameter regime as one of global chaos since the 
kinematically accessible phase space is predominantly chaotic. 

The dynamics of the classical and quantum distributions are much simpler in this regime. 
We find that initially localised distributions, launched from arbitrary initial conditions, re- 
lax to the microcanonical distribution on a very short time-scale. To demonstrate this, we 
consider the dynamics of an initial quantum state with s = 140 and l = 154, and a corre- 
sponding classical density, launched from 9{0) = (20°, 40°, 160°, 130°). Though the initial 
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Figure 4.4: Same as in the previous figure, but for PL^{mi) on the left and Pj^irrij) on 
the right, at time-step n = 200. Both Pj^rrij) and Pj^^rrij) are localised relative to the 
projected microcanonical distribution PJ^J^{mj). 

distributions are the same as in the mixed regime, by time-step n ~ 6 the quantum and 
classical distributions have aheady spread to the system dimension and begin to exhibit 
noticeable differences. By time-step n ~ 12 both distributions have relaxed very close to 
the microcanonical distributions. We plot the equihbrium quantum and classical projected 
distributions Phzifnj) and Pj^i'rnj) in Fig. for time-step n = 50. The projected classi- 
cal distributions are nearly indistinguishable from the microcanonical forms, PY^^{mj) and 
P'^^{mj), and the quantum distributions again exhibit small fluctuations about the classical 
distributions. We have found that these equilibrium quantum-classical differences asymptote 
to a non-vanishing minimum when the measure of KAM surfaces becomes negligible. These 
minimum quantum fluctuations reflect characteristic deviations from the microcanonical 
state that arise because the equilibrium quantum state is a sequence of pure states, whereas 
the microcanonical state corresponds to a random mixture. 
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Figure 4.5: The equilibrium shapes of PL^{mi) and Pj^{mj) at time-step n = 50 with 
l = 154 for a state launched in the global chaos regimc (7 = 2.835, r = 1.1, a = 5). The 
quantum distributions exhibit small rapidly oscillating fluctuations about the projccted 
microcanonical distributions. The classical distributions are not visible since the points lie 
within the fluctuating quantum data. 
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4.3 Rates of Relaxation to Equilibrium 

In order to characterize the time-scale of relaxation to equilibrium it is convenient to study 
the time-dependence of a scalar measure that is sensitive to deviations from the equilibrium 
state. A conventional indicator of this rate of approach to equilibrium is the coarse-grained 
entropy, 

H=-J2P^^^Pi■ (4.11) 

i 

Here the {Pi} stand for the quantum probabilities associated with projectors onto some 
basis of microstates {e.g. the projected distributions discussed in the previous section). The 
sum ( [4.11 ) is a Standard measure of the Information contained in a probability distribution 



and is sometimes called the Shannon entropy. 

The Shannon entropy has a number of useful properties. First, unlike the von Neumann 
entropy Tr[plnp], the Shannon entropy is basis-dependent. It reduces to the von Neumann 
entropy if the 'chosen' basis diagonalizes the state operator. However, this basis, or, more 
precisely, the set of projectors onto the (time-dependent) spectral decomposition of the state 
operator, does not necessarily correspond to a set of classically meaningful observables. Our 
main interest is to examine correspondence at the level of classical dynamical variables, so we 
consider probabilty distributions associated with projectors onto the eigenstates of classically 
well-defined operators. The classical counterparts to these probability distributions are 
associated with some fixed partioning of the phase space into celis of width h along the axes 
of the associated dynamical variable. 

Second, whereas the von Neuman entropy of the total system is constant in time 
(Tr[plnp] = since p = \^){ip\), the basis-dependent Shannon entropy may have time- 
dependence even if the quantum state is pure. Thus (4.11) may be applied to examine 



the rate of relaxation of either pure or mixed quantum states. It is in this sense that we 
use the term relaxation, although the time-evolution is unitary in the quantum model (and 
volume-preserving in the classical model). 

Given some fixed partioning of the phase space, if a classical state remains evenly spread 
through the phase space celis it occupies, and spreads through the phase-space exponentially 



with time, then an entropy like ( 4.11 ) should grow linearly with time. In this section we 
show that this argument holds approximately also for quantum states launched from a 
classically chaotic region of phase space. The actual rate of relaxation of the quantum 
states is accurately predicted by the classical entropy even for small quantum numbers. 
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Figure 4.6: Comparison of the quantum and classical entropies H[Jz] = 

~ (™i ) los ("^i ) for s = 140 and / = 154 in (a) regime of global chaos (7 = 2.835) ; 

(b) chaotic zone of the mixed regime (7 = 1.215); (c) regular zone of the mixed regime 
(7 = 1.215). 



We demonstrate this behaviour by first considering the quantum entropy Hq[Jz] of the 
probabilities associated with the eigenvalues mj of J^, i.e. the probabilities defined in (4.4). 
The corresponding classical entropy, Hc[Jz], is calculated using the discrete classical prob- 



abilities ( [4.81 ). In Fig. 4.6 we compare the time-development of the quantum and classical 
entropies using quantum numbers s = 140 and / = 154. For these quantum numbers, the 
microcanonical {i.e. maximum) value of the entropy is H^'^[Jz] = 6.2. In case (c), corre- 
sponding to a regular zone of the mixed regime (^(0) = (5°, 5°, 5°, 5°), 7 = 1.215), we actually 
see the greatest amount of difference between the quantum {Hq[Jz]) and classical {Hc[Jz]) 
entropies. Hg exhibits a quasi-periodic oscillation about its initial value whereas for Hc 
these oscillations eventually dampen. For smaller quantum numbers, and thus broader ini- 
tial states, Hc dampens much more rapidly although Hg continues to exhibit a pronounced 
quasi-periodic behaviour. In case (b), with initial centroid 9{0) = (20°, 40°, 160°, 130°) set 
in a chaotic region of the equally mixed regime, both Hg and Hc oscillate about an initially 
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Figure 4.7: Comparison of the quantum and classical subsystem entropies H[Lz] = 
— J2m PL^{n^i)^^PLA^i) increasing system sizes in the global chaos regime of Fig 



4.6 



increasing average before relaxing towards a constant value tl·iat lies well below the micro- 
canonical maximum H'^^[Jz] = 6.2. This saturation away from the maximum is expected 
in the classical model since a large fraction of the kinematic surface is covered with regular 
islands and remains inaccessible. In case (a), corresponding to the regime of global chaos 
(7 = 2.835) and with the same initial state as (b), the entropies are nearly identical. Both 
grow much more quickly than in the mixed regime case, roughly linearly, until saturating 
very near the maximum value. 

The quantum entropy is very well approximated by its classical counterpart also for 
smaller quantum numbers. In Fig. |4.7| we display the growth rates of the quantum and 
classical entropies of the probabilty distributions associated with the observable for three 
sizes of quantum system (/ = 11, / = 22, / = 220) using the same parameters and initial 



condition as for data-set (a) in Fig. 4^. In each case the quantum entropy is essentially 
identical to the corresponding classical entropy. The initial rate of growth is similar in each 
case, roughly linear, and of order the Lyapunov exponent, = 0.45. 

These results extend previous work demonstrating that the widths of quantum states 
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grow exponentially with time, on average, until saturation at the system dimension [46, 11, 



35 1. Modulo the small quantum fluctuations, for both quantum and classical models we find 



that the subsequent relaxation to an equilibrium configuration occurs on the time-scale, 

trel ~ tsat + ©(A^^), (4.12) 

where tgat — A^^^ln/ estimates the time it takes the initial coherent state to reach the 
system dimension. The exponent A^, is the exponent governing the growth rate of the state 
width |35|. The last term 0{X2^) approximates the additional time-required for the state 



to become more or less uniformly spread over the accessible phase space. In predominantly 
chaotic regimes we have found that A^ ~ A^ , though in mixed regimes is generally a few 
times larger than the largest Lyapunov exponent. 
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4.4 Time-dependence of Quantum- Classical DifFerences 

Before examing the scaling of quantum-classical difFerences with increasing quantum num- 
bers, it is useful to determine first their time-domain characteristics under the different 
types of classical behviour. Previous work has shown that quantum-classical differences for 
low-order moments, though initially small, grow exponentially with time when the classical 
motion is chaotic [pA| until the states approach the system size . On this saturation 
time-scale those quantum-classical differences reach their maximum magnitude, but surpris- 
ingly this maximum was small, 0{h). More specifically, it did not scale with the quantum 
numbers. Of course, two distributions can be altogether different even when the differences 
between their means and variances are quite small, and therefore it is useful to examine the 
differences between the quantum and classical states in a more sensitive way. 

In this section we examine the time-dependence of bin-wise deviations between the quan- 
tum and classical probability distributions. For the observable this indicator takes the 
form, 



1 ' 



^ [PLAm)-PiS^i)?- (4.13) 



\ (2/ + 1) — ; 

\ ^ ' mi=-l 

This Standard deviation estimates typical quantum-classical differences on the scale of h 
along the L^-axis. Each interval is centered on a quantum eigenvalue. The P^^{mi) cor- 
respond to a measurement, or coarse-graining, of the classical density on an extremely 
fine-scale. 



In Fig. 4.8 we examine the time dependence of (j[Lz] for the same three classical sets 
of parameters and initial conditions displayed in Fig. The initial value of (t[L^] is gen- 
erally not zero since it is not possible to match all the marginal distributions exactly in 
the case of the SU(2) coherent states The actual magnitude of the initial discrepancy 
depends on the angle between the axis of measurement, e.g. Lz, and the direction of po- 
larization of the initial state. For both chaotic states the differences initially decrease from 
their angle-dependent value and then increase until saturation at a steady-state value. This 
steady-state value is reached much later in the mixed regime (upper sòlid line), than in the 
global chaos regime (lower sòlid line). It occurs on the time-scale, trei, on which the un- 
derlying distributions have reached their steady-state configurations (modulo the quantum 
fluctuations) . 

As shown in Fig. |4.8|, the quantum-classical differences are actually largest for the regular 
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Figure 4.8: Time dependence of the Standard deviation cr[Lz] oi quantum-classical dif- 
ferences ( [4. 131) for states launched from a regular zone (dotted line) of the mixed regime 
(7 = 1.215), from a chaotic zone of the same mixed regime (middle sohd hne), and from 
the regime of global chaos (lower sòlid line, 7 = 2.835). The initial discrepancy is relatively 
large, but quickly decreases, and then increases until reaching an asymptotic equilibrium 
value. This occurs more slowly for the mixed regime case, for which the asymptotic value 
is also larger. In all cases s = 140 and / = 154. 



state (dotted line) of the mixed regime (7 = 1.215) at both early and late times (relative 
to the relaxation time-scale). The steady-state magnitude of the differences for the global 
chaos regime (7 = 2.835) is significantly smaller than the typical magnitude for the mixed 
regime. However, for larger vàlues of the classical perturbation strength 7, this average 
steady-state magnitude does not decrease further (with the quantum numbers held fixed) 
but has reached a non-vanishing minimum. The magnitude of the minimum steady-state 
fluctuations, ~ 2 x 10~^, should be compared with a typical magnitude of the quantum 
and classical distributions, Pi^(m;) ~ 3 x 10^^. In the following section we examine how 
these fluctuations scale with increasing quantum numbers. 

Above we have considered quantum-classical differences for observables (projectors onto 
subspaces) associated with the factor space Tíi. In this factor space the state is initially 
pure but becomes mixed as a result of dynamical interacions with the other subsystem. It 
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is interesting to check if the dynamical behaviours of the difFerences are an artefact of this 
dynamical mixing. Therefore we consider also bin-wise quantum-classical differences for an 
observable (J^ = Sz + L^) that acts non-trivially on the full Hilbert space Tí = Tís Tíi. 
The quantum state in the full Hilbert space remains pure throughout the time-evolution. 
We construct the same Standard deviation of the bin-wise differences between the quantum 
and classical probability distributions as above, 



1— It -r y^ïPj (mA - (m,)!^, (4.14) 

^ [2(s + /) + l)] -^'^ ^' •^-^ ' ^ ^ 



where mj G {/ + s, / + s — 1, . . . , — (/ + s)}. In Fig. O we compare a[Jz] in the same three 
classical regimes examined in Fig. [4.8| . Once again the regular state (dotted line) exhibits 
the largest quantum-classical differences, and the differences for both chaotic states (middle 
and lower sòlid line) grow to a steady-state value on the time-scale at which the underlying 
distributions relax to their equilibrium configurations. As above, the average value of the 
differences for 7 = 2.835 (lower sòlid line) correspond to a non-vanishing minimum, that is, 
the average value does not noticeably decrease for larger vàlues of 7. The minimum quantum 
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fiuctuations are again small when compared with the average height of the probability 
distribution, [2{s + l) + 1]"^ ~ 2 x 10~^. 

For 7 ~ 2.835, the measure of regular islands is already very close to zero and the 
classical system is nearly ergodic on V. Similarly, the quantum state is no longer constrained 
by any invariant classical structures but spreads almost evenly about the accessible Hilbert 
space. We find that the Standard deviations of the quantum fiuctuations that account 
for the equilibrium quantum-classical differences approach a non-vanishing minimum as 
the classical dynamics approach ergodicity on V. These equilibrium differences can not 
vanish (for fixed quantum numbers) because the total quantum state remains pure under 
the unitary dynamics, whereas the microcanonical equilibrium corresponds to an equal- 
weight mixture. 



CHAPTER 4. CORRESPONDENCE FOR THE PROBABILITY DISTRIBUTIONS 75 



4.5 Correspondence in the Classical Limit 

We now turn to an examination of the classical limit, J'/h — > oo, where is characteristic 
system action. Since the quantum-classical differences grow to their largest vàlues once the 
states have spread to the system size and subsequently fluctuate about this magnitude, we 
will examine the scaling of the differences in this late time-domain, that is, when the states 
have relaxed close to their equilibrum configurations. Moreover, these scaling results will 
then complement previous work that has focussed on correspondence at early times [^] , in 
the Ehrenfest regime when the states are narrow relative to the system dimensions. 

We wish to determine if the Standard deviation of the quantum-classical differences 
(defined in the previous section) decreases in magnitude with increasing quantum numbers. 
When comparing models with increasing quantum numbers, we hold the width of each 
probability bin fixed (at h = 1). Since the number of bins will increase with the quantum 
number, it follows that the height of the probability distribution in a given bin will also 
decrease. Consequently, we construct a scale-independent, or relative, measure of the bin- 
wise quantum fluctuations by taking the ratio of the Standard deviation to the average value 
of the probability distribution. For the observable this takes the form, 

R[L,{n)] = ^Í^M = Ni a[L,{n)]. (4.15) 

where the average value Pl^ = 1/(2Z-|-1) = 1/Ni. If this relative measure approaches zero in 
the classical limit then the quantum probabilty distribution converges to the corresponding 
classical one in that limit. 



In Fig. 4.10 we consider typical equilibrium vàlues of R[Lz{n)] plotted against l/^/Wi. 
We study the scaling using Ni because it is equal to the dimension of the factor space 
Tíi and it is also proportional to the subsystem size A'^; ~ 2|L|. We first consider a state 
launched in the global chaos regime (7 = 2.835, r = 1.1), with initial condition 9{Q) = 
(45°, 70°, 135", 70°). The scatter of plus signs for each A'^; = 2/ + 1 value corresponds to 
time-steps n such that 41 < n < 50. These time-step vàlues are chosen because they occur 
well after the relaxation time í^e/ — 6. In this regime the data exhibits very little scatter. 
A least-squares fit to the curve R = A/^/Nl + B yields a value for the intercept B that is 
consistent with zero {B = 0.001 ± 0.001) and a slope of order unity {A = 1.032 ± 0.02). 
An intercept consistent with zero implies that quantum-classical differences vanish in the 
classical limit, i.e. PL^{mi) — > Pj^^{mi) as / 00. This result is especially remarkable since 
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Figure 4.10: Scaling of relative quantum-classical differences (4.15 in the equilibrum time- 
domain versus increasing system size. Scatter of crosses corresponds to time-steps 191 < 
n < 200, for a state launched in the chaotic zone of the mixed regime (7 = 1.215). Scatter 
of plus signs corresponds to time-steps 41 < n < 50, for a state launched in the global chaos 
regime. Data sets in both of these regimes are consistent with the scaling law R ~ Ni~^^'^, 
where Nj = 21 + l. 



we have considered the differences that arise given classical measurements which resolve the 
observable Lz with the rather extraordinary precision of ?i = 1. 

We next consider a state launched from the chaotic zone of the mixed regime (7 = 1.215, 



r = 1.1) with é'(O) = (20°, 40", 160°, 130°). The scatter of crosses in Fig. |l0| corresponds to 
time-steps 191 < n < 200 n, again chosen well after the relaxation time trei for the range of 
quantum numbers considered. The scatter of quantum-classical differences at each A'; value 
is much more significant in this regime in which the equilibrium distributions reflect a much 
more complex phase space structure. However, the relative differences exhibit, on average, 
a similar dependence on the quantum numbers as in the predominantly chaotic regime. 
In this regime a least-squares fit to the curve R = A/ Ni^^"^ + B yields a slope of order 
unity {A = 3.39 it 0.15) but a negative value for the intercept {B = —0.017 it 0.009) within 
two-standard deviations of zero. A negative intercept is not physically meaningful (since R 
is a positive definite quantity) and we assume it arises as a consequence of the statistical 



CHAPTER 4. CORRESPONDENCE FOR THE PROBABILITY DISTRIBUTIONS 77 




1/2 

scatter in the data. Also plotted is the curve R = C/N^ , with slope C = 3.09 it 0.04 also 
determined from a least-squares fit. Both fits are good, with reduced vàlues of order 
unity. 

As we noted in the last section, the subsystem states do not remain pure, because of 
dynamically induced entanglement between the subsystems. Since the subsystem state ( [4.3D 
in the factor space Tíi is not pure, but highly mixed in the equihbrium time-domain, it is 
possible that the scaling with Ni that we observe is related to the purity-loss from this 
entanglement. Consequently, it is useful to examine the scaling of the quantum-classical 
differences for the total spin J^. The operator acts non-trivially in the full Hilbert space 
Tí. In this Hilbert space the system is described by a pure state vector at all times. In Fig. 



4.11 we consider the scaling of the ratio, 

R[Un)] = = N, a[Un)], (4.16) 

where Pj, = [2{s + /) + 1]"^ = iVri is the average value of either distribution, versus the 
dimension Nj. Here Nj is the number of subspaces associated with distinct eigenvalues (rrij) 
of the quantum operator (Jz)- In contrast with A'^;, Nj is not equal to the dimension of the 
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corresponding Hilbert space, though it is a measure of the system size since A',- ~ 2|J|. The 



parameters and initial conditions shown in Fig. 4.11 are the same as in Fig. 4.1C. The same 

1/2 

fit procedure as above, but appUed to the function R = A/N- + i?, yields a value for B 
that is again consistent with zero {B = 0.00038 ±0.0016) and a positive slope of order unity 
{A = 2.00 ib 0.04) in the predominantly chaotic regime (scatter of plus signs). Thus the 
relative Standard deviation for also decreases as the square of the quantum numbers and 
fits to an intercept that is consistent with zero. This imphes that the fluctuating quantum 
distributions approach the classical equilibrium, even for a few degree-of-freedom system, 
which is described at all times by a pure state. In a chaotic state of the mixed regime 
(scatter of crosses), the fluctuations are larger, and the same fit procedure as above gives 
[B = — 0.016ib0.012, A = 6.4ib0.3), where the negative value for B lies within two Standard 
deviations of zero and is presumed to result from the statistical scatter of the data. Also 

1/2 

plotted is the equation R = C /N- , with C = 5.97 it 0.06 determined from a least-squares 
fit. The fits to both equations are good, with reduced vàlues of order unity. 
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4.6 Discussion 

We have shown that, in classically chaotic regimes, initially localised quantum states relax to 
an equilibrium configuration that reflects the details of the classical phase-space structure. 
We find a remarkable degree of correspondence between the quantum and classical relaxation 
rates, even for small quantum numbers. Moreover, contrary to results obtained for the low- 
order moments |jïï|, the degree of difference between the probability distributions is 
actually smaller for the chaotic states than the regular states. 

The equilibrium quantum distributions exhibit small rapidly oscillating fluctuations 
about the coarse-grained classical equilibrium. As the measure of regular islands on the 
classical manifold approaches zero, the quantum and classical equilibrium configurations 
approach their microcanonical forms, and the quantum fluctuations about the classical 
equilibrum approach a non-vanishing minimum. This minimum arises because we consider 
total quantum states that are pure, whereas the microcanonical configuration is produced 
by an equal-weight mixture. 

For the distributions associated with the subsystem observable L, the Standard deviation 

— 1/2 

of these differences, relative to the average value, decreases as , where A'^; = 2/ + 1 ~ 

2|L| is the dimension of the factor space, and becomes vanishingly small in the limit of 
large quantum numbers (i.e. large spins). These results suggest that correspondence with 
classical Liouville mechanics emerges in the classical limit for time-scales much longer than 
the Ehrenfest time. 

A great deal of recent work has emphasized that the loss of purity resulting from inter- 
actions with a quantum environment removes characteristic quantum effects and improves 



the degree of quantum-classical correspondence |58, 56, 52 1. While this is certainly the 
case for small quantum systems, it has been further argued that these decoherence effects 
must be taken into consideration to see the emergence of classical properties from quantum 
mechanics, even in the limit of large quantum numbers, if the classical motion is chaotic 

Since our model is comprised of interacting subsystems, initially separable pure states 
become entangled dynamically; the subsystem states (in each factor space) do not remain 
pure but become mixed. This entanglement process has an effect that is analogous to the 
process of decoherence. Hence one might suspect that the emergent classical behaviour 
that we have observed for the properties of the subsystem L may be strictly the result of 
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a "decoherence" efFect arising from entanglement with the other subsystem. To address 
this possibility, we have considered also the quantum-classical differences that arise in the 
probability distributions for a total system observable, J^. In this case the quantum observ- 
ables are projectors onto subspaces of the full Hilbert space, rather than merely a factor 
space. The quantum state in this full Hilbert space is not subject to any entanglement 
or decoherence and remains pure throughout the unitary time-evolution. We have found 
that the scale-independent Standard deviations for these quantum-classical differences de- 
crease as l/^/77j, where Nj = 2{s + l) + 1 ~ 2|J| is a measure of the system size and 
N = {2s + 1)(2Z + 1) is the dimension of the Hilbert space. The bin-wise quantum-classical 
differences become increasingly difficult to observe, in the limit of large quantum numbers, 
even for system observables that are isolated from the effects of decoherence. In this sense 
the process of decoherence is not necessary to produce quantum-classical correspondence in 
the classical limit. 



Chapter 5 

Breakdown of Ehrenfest 
Correspondence 

5.1 Introduction 

The quantum-classical correspondence principle expresses the view that the predictions of 
classical mechanics must emerge from quantum mechanics in the macroscopic hmit. A typ- 
ical statement of this principle is given by Messiah [^]: "In the hmit ?i ^ 0, the laws 
of Quantum Mechanics must reduce to those of Classical Mechanics." Here the expres- 
sion "h 0" is shorthand notation for the limit where the characteristic system action is 
much larger than Planck's constant. Quantum-classical correspondence is required in the 
macroscopic limit if classical mechanics provides a vàlid approximation to the observed phe- 
nomena in this limit, whereas quantum mechanics provides correct predictions on all action 
scales. In order to test whether this principle holds in relevant physical situations, it is 
necessary to make more precise statements about which classical theory of mechanics, or, 
put differently, which classical properties, should be expected to emerge from the underlying 
quantum description. In this context, it is useful to introduce and distinguish two different 
interpretations of the quantum-classical correspondence principle. 

The first of these, which I shall call Ehrenfest correspondence, is the claim that the 
phase space trajectòries described by Hamilton's equations of motion are required to emerge 
from quantum mechanics in the macroscopic limit. This approach to characterizing corre- 
spondence was originally devised by Ehrenfest |^] and forms the basis for most textbook 
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discussions of the correspondence principle |65, 81, 68, 88, 67]. For a genèric Hamiltonian 
system of the form H = /2m + V{q), this theorem states that the time-dependence of 
the expectation vàlues {q{t)) = {il){t)\q\ilj{t)) and {p{t)) = (V'(í)lí'IV'(í))) in an arbitrary 
quantum state '0(í)> is prescribed by the differential equations, 

dm) ^ (pM 

dt m 

= (Fh». (5.1) 

If the quantum state is sufíiciently well-locahsed, then the approximation, 

{F{q)) F{{q)) (5.2) 

holds, and the differential equations for the quantum expectation vàlues are then well ap- 



proximated by Hamilton's equations of motion [51 1. A clear statement of this interpretation 
of the correspondence principle is articulated by Merzbacher [^]: "We require that the clas- 
sical motion of a particle be approximated by the average behaviour of a wave packet with 
a fairly sharp peak and as precise a momentum as the uncertainty principle permits and 
that the expectation vàlues of the dynamical variables, calculated for such a wave packet, 
satisfy the laws of classical mechanics." 

The second interpretation of the correspondence principle, which 1 will refer to as Liou- 
ville correspondence, is the weaker proposition that only the statistical properties of Liouville 
mechanics are required to emerge from the quantum description in the macroscopic limit. 
This view may be characterized by the condition that the quantum expectation vàlues, 
{A{q,p)), for classically well-defined observables, should approach classical ensemble aver- 
ages for the associated dynamical variables, {A{q,p))c. In this correspondence picture, the 
classical ensemble averages are calculated from the prescription, 

{A)c = JdqJ dpAiq,p)p,{q,p,t) (5.3) 

where the classical density, Pc{q,p), describes the possible configurations for a system with 
incompletely specifed phase space coordinates. This density evolves with time according to 
the Liouville equation, 

dp{q,p,t) 

g-^ = {p[q,P, t),H} (5.4) 

where {•, •} is the Poission bracket. In the previous chapters I have examined the character- 
istics of Liouville correspondence in considerable detail. I have demonstrated, for the model 
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system described in Chapter 2, that the statistical properties of Liouville mechanics emerge 
from quantum mechanics with a degree of approximation that improves as the system action 
increases. In particular, I examined this correspondence for quantum expectation vàlues and 
classical ensemble averages associated with the system dynamical variables, their variances, 
and characteristic functions defined over extremely small intervals along the phase space 
axes corresponding to these variables. Even for the chaotic states of the model, the differ- 
ences between quantum and Liouville mechanics for these observables have been shown to 
become vanishingly small in the limit "?i — > 0" over physically relevant time-scales. Most 
significantly, the correspondence with Liouville mechanics remained vàlid not just at early 
times (for well-localised quantum states), but also at late time (well after the width of the 
quantum state had grown to the system size). The conditions of Liouville correspondence 
are therefore satisfied in the macroscopic limit, in the sense that the observables of quan- 
tum mechanics are well approximated by the statistical predictions of classical Liouville 
mechanics, over physically accessible time-scales, even in the case of chaotic motion. 

In this chapter I will demonstrate that the conditions of Ehrenfest correspondence may 
not be satisfied in the macroscopic limit when the classical motion is chaotic. In particular, I 
will argue that Ehrenfest correspondence may fail on physically relevant time-scales for some 
chaotic macroscopic systems, and, therefore, that quantum mechanics may be unable to 
describe the observable, deterministic motion of some macroscopic bodies. The implications 
of this result for the interpretation of the quantum state will be drawn out in the discussion 
at the end of this chapter. 

This chaper is organized as follows. First I will present a theoretical argument indicating 
that the presence of chaos in general Hamiltonian models leads to macroscopic differences 
on a time-scale that grows only logarithmically with increasing system action. Specifically, 
I will show the quantum state centroids deviate from the predicted Newtonian trajectory 
in an experimentally observable manner. This time-scale may be experimentally accessi- 
ble for some real macroscopic systems, leading to an observable breakdown of Ehrenfest 
correspondence. I will then illustrate the features of this general argument in the specific 
case of the coupled spin system described in Chapter 2. This example will also clarify how, 
contrary to a recent claim in the literature [^], the effects of decoherence are unable to 
restore this anticipated breakdown of Ehrenfest correspondence. In the discusssion I will 
explain how, on the assumption that individual macroscopic systems conserve energy and 
remain well-localised over experimentally verifiable time-scales, this argument leads to the 
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conclusion that quantum mechanics is unable to provide a vàlid description of macroscopic 
dynamics, and is therefore not a complete theory. 

5.2 Ehrenfest Correspondence Conditions 



The conditions for Ehrenfest correspondence are summarised succinctly by Messiah |67]: 
"In order that this picture be satisfactory, it is necessary that: (a) the mean vàlues follow 
the classical laws of motion to a good approximation; (b) the dimension of the wave packet 
be small with respect to the characteristic dimensions of the problem, and that they remain 
small in the course of time." It is useful to formulate these requirements as mathematical 
conditions. 

Let Xc = {ql, . . . ,p},, . . . , ) , denote the 2A'' phase space coordinates of an N degree 
of freedom system, let (x) = ('í/'(í)|x|'0(í)) denote the corresponding quantum expectation 
vàlues, and let = ((x*)— x*) stand for the Ehrenfest differences. Further, let denote some 
prescribed set of thresholds that characterize adequate agreement between the quantum and 
classical predictions (e.g., determined by the resolution of the experiment measurements) , 
and let toj,s stand for a characteristic time-scale over which the system may be subject to 
experimental observation. I will use to denote a characteristic action of a physical system, 
with the macroscopic limit characterized by a very large value of the system action relative 
to Planck's constant, J jfi » \. In order for Newtonian trajectòries to emerge from the 
quantum description, the Ehrenfest differences must remain small relative to the prescribed 
threshold, 

\Ç{t)\ = \(x^) - xl\ < x" for t < tobs and J /h » 1. (5.5) 

For a genèric Hamiltonian system, this is only possible if the quantum state remains sufíi- 
ciently well-localised, 

Ax* < X* for t < tobs and J/h » 1, (5.6) 

where for notational convenience I have used the same tolerance thresholds for the quan- 
tum state widths Ax^ = [((x*)^) — (x*)^]^/^ as for the Ehrenfest differences. These Ehrenfest 
conditions may be taken to define the proposition that quantum mechanics can describe the 
deterministic motion of a single macroscopic system, to within some prescribed accuracy and 
over physically observable time-scales, through the centroids of a well-localised state. I will 
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now show that these Ehrenfest conditions may break down on physically relevant time-scales 
for macroscopic bodies that are subject to classically chaotic dynamics. 

5.3 The Ehrenfest Break-Time 

Let the Ehrenfest break-time, tEhr, be defined as the shortest time upon which one of the 



Ehrenfest differences grows larger than one of the thresholds in Eq. (5.5). The question 
under consideration is whether tEhr < tobs holds in relevant macroscopic situations. 

A simple argument sufíices to estimate how the time-scale ÍEhr grows with with in- 
creasing system action J'. In the general case of an autonomous classical flow, the time- 
dependence of the 2N phase space variables are prescibed by the following set of first-order 
differential equations, 

r/x 

^ = {x„/í}=F(x,). (5.7) 

where {•, •} is the Poisson bracket and H stands for the Hamiltonian. The time-evolution of 
the corresponding Heisenberg operators, Xg, is given by the same set of differential equations, 



dxq i 



' - ^[x„i/]=F(x,), (5.8) 



where [•, •] denotes the commutator for two operators. Expanding the function F{xq) about 
the classical trajectory, Xc, gives, 

í7x f)V 

= F(x,) + (4 + 0{[xÍ - xif). (5.9) 

Taking the expectation value of both sides it follows that, 

While the quantum-classical differences remain sufficiently small, 

^ ^ dxjdxk dx'· ^ ' ^ ^ 



the higher-order terms in the expansion ( |5.9D remain negligible, and the growth of the 
Ehrenfest differences is governed by the matrix. 
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which is just the classical tangent map. The time-dependent eigenvalues of this matrix, 
evaluated along a fiducial trajectory, determine the stabihty of the classical flow. Conse- 
quently, in the case of classically regular motion, the quantum-classical differences will grow 
at most as a small power of time, 

m ^ Cot", (5.13) 

where a ~ 1. On the other hand, in classically chaotic regions, the local flow will exhibit 
expansion along some directions and contraction along others, and the very definition of 
the classical Lyapunov exponents |Q entails that the quantum-classical differences should 
grow, on average, exponentially with time, 

Ç(í) ~Çoexp(ALÍ), (5.14) 

where is the largest Lyapunov exponent. I have used the symbol Ço to denote the non- 
vanishing difference that will immediately arise between the quantum expectation vàlues 
and classical dynamical variables due to the presence of a non-zero quantum-state width in 
either position or momentum Q. Since denotes the threshold that defines a break between 
the quantum and classical predictions, it follows that this break arises on the time-scale, 

tEhr-^Xï'HxyO- (5-15) 

The result ( |5l5| ) can be expressed in a slightly more useful form. The "classical limit" 
of quantum mechanics corresponds to a sequence of quantum models that describe physical 
systems with increasing size (e.g., increasing quantum numbers), but with other parameters 
adjusted so that each model of this sequence is associated with the same classical system. 
The duration of the correspondence between the quantum and classical expectation vàlues 
is optimized if the initial states are always chosen to be minimum uncertainty states (i.e., the 
initial variances are not increased in proportion to the system size but held fixed), whereas 
the break thresholds are always chosen as fixed fraction of the characteristic system size 
for each model in this sequence. Under these conditions, the ratio (x^/C) ^^^^ scale as a 
power of the ratio (J' /h), where ^ is a characteristic system action, and the break between 
the quantum and classical theories arises on the time-scale, 

tEhr-- Xl^HJ/^), (5.16) 

which estimates the optimal duration of Ehrenfest correspondence for chaotic systems. The 
time-scale (5.16) has been derived previously for a number of specific model systems |15, 16, 
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17, 95 



54] through a variety of diíFerent techniques. Although some of these authors 



have interpreted the break-time ( ^.16| ) as a time-scale during which "classical" behaviour 
emerges from the quantum description 25, 48, ^4], 22 1, this interpretation is imprecise 
[0, 35, 36|, and I will refer to ( 5.16D specifically as the Ehrenf est hieak-time since it is defined 
as the time-scale during which the Ehrenfest correspondence conditions (|5.5|) remain vàlid. 



5.4 The Breakdown of Ehrenfest Correspondence 

Some authors have maintained that the break-time scahng in Eq. (5.16) is compatible with 
the requirement of quantum-classical correspondence since ÍEhr — > oo as "7i — > 0" |95, 



25, 22|. Of course, for real physical systems h is a constant, and the relevant question is 
whether (5.16) is sufíiciently long to accomodate the Ehrenfest correspondence conditions 
for all chaotic macroscopic systems. Assuming the thresholds designate macroscopically 
observable differences, then the ratio (xVO '^^ course, enormous. However, Zurek and 



Paz |9^, 99| have noted that the factor \og(J'/h), appearing in the break-time expression 
( 5.16| ), may be quite small for some macroscopic systems, even ones with astronomically large 
vàlues of the ratio (J^/h), since this ratio appears as an argument inside the logarithm. 
Consequently, the Ehrenfest break-time may be short compared to actual time-scales of 
observation. As an explicit example, these authors consider Hyperion, a moon of Saturn 
which is believed to exhibit a chaotic tumble [94, BO]. They estimate the upper bound í ~ 20 



years for the onset of a massive discrepency between the quantum and classical predictions. 

Zurek |9£] has suggested that such a short break-time raises the possibility of an observ- 
able breakdown of the correspondence principle for classically chaotic systems. Although 
I will demonstrate below that macroscopic chaotic motion certainly raises the possibility 
of an observable breakdown of the Ehrenfest correspondence conditions, it is important to 
emphasize that these situations do not suggest a breakdown of the correspondence princi- 
ple, since the possibility of observing macroscopic differences between quantum mechanics 



and classical Liouville mechanics appears to be extremely remote |35, 36 1. In the following 
section 1 will draw out some of the implicit assumptions of this breakdown argument that 
will help clarify this distinction. 
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5.5 Macroscopic Deviations and Kinematic Constraints 



The argument leading to ( 5.141 ) leaves room for significant deviations from the exponential 



growth rate over short times, since the Lyapunov exponent is strictly defined as an average 
taken in the Umit t — > oo. However, the exponential growth rate ( |5.14| ) and the log{J'/h) 
break-time should hold with a degree of approximation that improves as the initial differ- 
ence Ço decreases relative to the characteristic system size. A more critical assumption of 
the breakdown argument is that Eq. ( |5.16| ) describes the onset of macroscopic Ehrenfest 



differences. This follows only if the exponential growth rule ( 5.14 ) remains vàlid until the 
Ehrenfest differences grow to a significant fraction of the system size. But it is clear that 
the exponential growth of the Ehrenfest differences must eventually become invalidated. 



For example, in the case of bounded systems, Eq. (5.14) can remain vàlid at most until 
the differences grow to the size of the accessible phase space, at which time the growth 
will completely saturate. Although from the condition ( ^.11 ) it appears that any charac- 



teristic threshold at which this growth rate terminates may be expected to scale with the 
system size, in the absence of a rigorous argument, it is necessary to confirm this expectation 
through explícit calculation in specific model systems. Moreover, for any given system, it 
is necessary to determine if the Ehrenfest differences grow to a sufficiently large fraction of 
the system size that these differences become larger than the resolution of the macroscopic 
measurements. 

It may be objected that, even if the Ehrenfest differences grow to a macroscopic size 
on a physically relevant time-scale, these differences may not be experimentally significant 
as a result of the extreme sensitivity to initial conditions exhibited by the chaotic classical 
dynamics. Let Ag(0)Ap(0) ^ denote the non-vanishing area of phase space from which 
the classical trajectory is known to have originated. From a purely classical point of view, 
at some later time t, it will be impossible to predict the location of the classical trajectory 
to within an area smaller than Aq{t)Ap{t) ~ Ag(0)Ap(0) exp(2ALÍ) as a result of the 
exponential growth of the imprecision in the initial phase space coordinates. It is therefore 
impossible to experimentally confirm the Newtonian prediction with a precision better than 
Ag(0)Ap(0) exp(2ALÍ). However, according to (|5.14 ), the Ehrenfest differences are expected 



to grow at the exponential rate, ÇgÇp ~ Çq(0)Çp(0) exp(2ALÍ), which is the same rate of 
growth governing the imprecision of the classical prediction. Therefore the correctness of 
the classical (Newtonian) theory may not be experimentally confirmed with a precision that 
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is smaller than the Ehrenfest differences. According to this objection, even if the predictions 
of the quantum and classical theories are macroscopically difFerent, these differences are 
experimentally indistiguisl·iable, and tfie tfieoretical prediction of a macroscopic discrepancy 
does not indicate an observable breakdown of the Ehrenfest correspondence principle. 

However, this practical objection may be overcome by considering certain kinematic 
constraints satisfied by the classical dynamical variables. For an autonomous system with 
energy £^(q, p), although a prediction of the classical dynamical variables q(í) and p(í) 
remains subject to an exponentially growing imprecision, these variables must satisfy the 
time-independent constraint £^(q, p) = Eg, since the energy is a constant of the motion. 
Consequently, at any given time, the experimentally measured vàlues of the coordinates, 
q,Tn{t) i and Pm(í) =t Sp, are predicted to satisfy the constraint £'(q^(í), p^(í)) = Eq on 
the basis of the classical theory. 

On the basis of the quantum theory, the expectation value of the energy (-^(q, p)) = Eo 
is also a constant of the motion since the energy operator commutes with the Hamilto- 
nian. However, the quantum mechanical prediction for the time-dependence of the function 
£'((q(í)), (p(í))) is not subject to the same constraint as the classical function, and may ex- 
hibit macroscopic deviations from that constraint. Consequently, if one can show that the 
quantum prediction E'((q(í)), (p(í))) differs from the classical prediction £'(q(í), p(í)) = Ef, 
by a macroscopically large magnitude, then experimental observation of the phase space 
coordinates (with adequate precision) will distinguish the predictions of quantum theory 
from thosc of the classical (Newtonian) theory. If the function -E(q, p) is reasonably well- 
behaved, then a set of macroscopically large Ehrenfest differences, arising on the time-scale 
tEhr, should lead to a macroscopically large deviation between the quantum prediction, 
E{{ci{tEhr)),{p{tEhr))), and the classical one, E{cí{tEhr),P{tEhr)) = K- This argument 
may be generalized to the case of a non-autonomous system provided the system possesses 
any constant of the motion, /(q, p), i.e., a function that satisfies. 



dljcí, p) 

dt 



= {-^(q, p),-f^(q,p, 



t)} = 0. 



(5.17) 



5.6 Case Study: Coupled Spins 



I will examine these features of Ehrenfest correspondence by explícit calculation of the 
quantum and classical dynamics for the model of nonlinearly coupled spins described in 
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Chapter 2, 

i/ = a(5, + L,)+c5,L,^ó(í-rn). (5.18) 

n 

The key feature of the theoretical argument of section that must be checked by explícit 
calculation is whether the exponential growth of the Ehrenfest differences persists until 
reaching a magnitude that scales with the system size. 

I will consider this correspondence in the predominantly chaotic regime associated with 
the classical parameters 7 = 2.835, r = 1.1, and a = 5, and with largest classical Lyapunov 
exponent = 0.45. The quantum dynamics are calculated using an initial coherent state 
centered at 0(0) = (20°, 40°, 160°, 130°). The initial Ehrenfest differences may be set to zero 
if the classical initial conditions are set equal to the initial quantum state centroids, 

|L(0)| = |(L(0))|=/. (5.19) 

In the case of the coordinate L^, the time-dependence of the Ehrenfest difference, 

CL,it) = {Uit))-Li{t), (5.20) 



is plotted in Fig. 5.1 for quantum numbers / = 22 and / = 220 on a semi-log scale. The 
upper envelope of the Ehrenfest differences indeed corresponds to approximately exponential 
growth, as expected from Eq. ( 5.14| ). More importantly, this exponential growth persists 
until saturation at the system size that is associated with each quantum model. The system 
sizes corresponding to / = 22 and / = 220 are indicated by the horizontal sòlid lines in 
the figure. These results are consistent with the general picture of the chaotic quantum 
dynamics developed in Chapters 3 and 4, where it was shown that the expectation vàlues 
of dynamical variables approach their microcanonical equilibrium vàlues, e.g., {Lz) — > 0, on 
the saturation time-scale. 



As explained in section even macroscopic vàlues of Ehrenfest differences in the co- 
ordinates, such as ( [5. 201 ), arising on a physically accessible time-scale, do not automatically 
imply the possibility of an observable (experimentally significant) breakdown of correspon- 
dence, due to the exponential loss of precision that characterizes the chaotic classical dy- 
namics. This objection may be overcome by considering the constraints imposed on the 
dynamics by the classical constant of the motion. For the model ( ^.18 ), the magnitudes of 
the subsystem spins provide the required constants of the motion. 



CHAPTER 5. BREAKDOWN OF EHRENFEST CORRESPONDENCE 



91 



103 F 



102 r 



{L'z — {Lz 



0.1 



10 r- 



# • • • • 



"o cr 

o 



o 



s=20, 1=22 o 
s=200, 1=220 • 



10 



12 



14 



Kick Number 



Figure 5.1: The time-development of the Ehrefcst diffcrcnccs for Lz for classical parameters 
7 = 2.835, r = 1.1 and o = 5 for which the phase space is predominantly chaotic. Open 
circles correspond to quantum numbers (s, l) = (20, 22) and filled circle correspond to (s, /) = 
(200, 220). In both cases the approximately exponential growth of the Ehrenfest differences 
persists until saturation at the system size |L|, indicated by horizontal hnes. 



^ = {S^H}=0. (5.21) 

As a result of thesc constant functions, the time-dependent spin components are subject to 
the time-independent contraints. 

Llit) + Llit) + Llit) = |Lp 

Slit) + S'yit) + Slit) = |Sp (5.22) 

In addition, the magnitudc of the total system angular momentum, though time-dependent, 
is subject to the constraint, 

L2 + s2 - 2|L||S| < J2 < + + 2|L||S|. (5.23) 

In the case of the quantum model, the cxpcctation vàlues of the operators and 
are conserved, but the associatcd quadratic functions, 

(L)2 = (L,(í))2 + (L,(í))2 + (L,(Í))^ 
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Figure 5.2: The time-development of the differences between the classical constant of the 
motion = P and the quantum prediction for (L)^ for the same classical parameters as in 
Fig. 5.1. The data are compared to an exponential growth rate given by twice the classical 
Lyapunov exponent, A^^ = 0.45. Open circles correspond to quantum numbers (s,/) = 
(20,22) and filled circle correspond to {s,l) = (200,220). In both cases the exponential 
growth of the Ehrenfest differences persists until saturation at the maximum possible value 
L^, indicated by horizontal lines. 



{Sf = {S,(t)f + {Sy{t)f + {S.it)f, (5.24) 
are not conserved quantities. Similary, the magnitude of the total angular momentum, 

(J)2 = (J^)2 + (J^)2 + ^J^)2^ (5^25) 



is not bounded by the classical constraint ( 5.23| ). 

In Fig. ^ the time-dependent differences L2 - (L)2 are plotted for í = 22 and í = 220. 
The growth of these differences is approximately exponential, and governed by a rate that 
appears to be independent of the system size. This growth may be estimated on the basis of 
the analytical argument predicting exponential growth ( |5.14[ ) for the Ehrenfest differences 
for the coordinates, Ç = L — (L). Writing L(n) = (L(n)) + ^(n) and squaring both sides 
gives, 

L'^ - (Lf = f - 2^- (L). (5.26) 
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Since ~ exp(2ALn) grows very rapidly, and assuming the second term in (5.26) may 

therefore be neglected for large n, it follows that 

L2_(L)2^Ç2g^p(2ALn). (5.27) 

An alternative approach to characterizing the growth rate for the Ehrenfest difference 
— (L)2 follows from comparison with the numerically measured exponential rate of growth 
of the quantum variance ( p.!] ). The magnitude of quantum variance (AL)^ provides a good 
approximation to the time-dependent Ehrenfest difference, — (L)^, for large l, that is, 

(AL)2 ~ - (L)2, (5.28) 

since, 

= /2 ~ (l2) = ;(; + l), (5.29) 
where quantities on the left correspond to the classical magnitude of the spin and those on 



the right to the quantum magnitude. Therefore, from ( 5.2g| ) and ( |3.1[ ), it may be expected 



that the Ehrenfest differences for the quadratic functions grow according to, 

- (Lf ~ / exp(2A^n), (5.30) 
Ai in regimes of predominantly chaotic dynamics (as shown in Chapter 3). This 



estimate of the growth rate is the same as (5.27) obtained on the basis of the analytical 



argument leading to ( |5.14 ), but also provides an estimate of the prefactor ~ L The 
exponential rate (|5.30| ) is plotted in Fig. for / = 22 and / = 220, with Xw = = 0.45, 
which slightly underestimates the growth rate of the data. 

The most important feature of the data in Fig. |5.2| is that this exponential growth of these 
observable Ehrenfest differences does not saturate until the differences reach the magnitude 
of the corresponding classical constant of the motion L^. The classical magnitudes = 484 



and ~ 48400, for L = 22 and L = 220, respectively, are plotted as sòlid lines in Fig. 5.2 . 
These are the maximum vàlues that may be expected on purely kinematic grounds, since, 

< |(L)2 -L^l < L^. (5.31) 

The break-time for the experimentally observable differences — (L)^ is determined 
numerically from the earliest time at which these differences exceed a fraction / of the 
classical magnitude, i.e., the condition, 

- (L)2| > /L^. (5.32) 
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Figure 5.3: The break time obtained by the condition ( |5.32| ), using a threshold /L^ = 
0.25L^, plotted against quantum number / ranging from 11 to 220. The data are obtained us- 
ing the same classical parameters as in Fig. 5.2 and initial condition 6 = 20°, 40°, 160°, 130°). 
The sohd hne corresponds to Eq. 5.33 with exponent A = 0.51 obtained by a least-squares 
fit. The dotted hne corresponds to Eq. 5.33 using the largest Lyapunov exponent Xl = 0.45. 



This calculated break-time is plotted against increasing system size, using / = 0.25, in Fig. 
( |5.3D . The numerical data exhibit a step-wise behaviour since the time, n, is coníined to 
integer vàlues. This scaling behaviour of the break time may be estimated by setting the 
LHS of the analytical prediction ( ^.301 ) equal to the threshold /L-^, and solving for the time, 
to give, 

tEhr^l/{2X^)log{fL). (5.33) 



This curve is plotted in Fig. 5^ using Xw = Xl = 0.45 (dotted line), which slightly over- 
estimates the break-time results. A least squares fit to ( p.33D gives Xw = 0.51, which 
provides an excel·lent fit to the data (sòlid line). The scaling demonstrated by the break- 
time data in this figure corresponds to that of a typical initial coherent state evolved under 
parameters that produce a predominantly chaotic classical phase space. These results ex- 
plicitly confirm that the log{J /h) dependence of the Ehrenfest break-time holds even for 
experimentally observable differences that have grown to a significant fraction of the relevant 
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classical magnitude. 



5.7 Discussion 

The argument outlined above should apply to a wide variety of chaotic systems. In particu- 
lar, the quantum expectation vàlues may be experimentally distinguished from the classical 
coordinates on the short time-scale ( [5.16| ) provided the classical model admits a constant 
of the motion which is some "reasonable" function of the coordinates. For example, the 
differences between the quantum and classical predictions for any polynomial function of 
the coordinates should grow exponentially as long as the differences between the coordinates 
are themselves subject to exponential growth. 

Though I have considered this problem in a specific model system, there is also evidence 
in other model systems that the break-time ( 5.16| ) should apply to Ehrenfest differences of 



order the system size when the phase space is predominantly chaotic. Fox and coworkers 



[47, 46 1 have shown that the quantum variances in the kicked top and the kicked rotor 
grow at an approximately exponential rate which subsides only when the variances have 
reached the system size. By virtue of the connection between the quantum variance and 



the time-development of the expectation vàlues expressed in (1.1), it follows that these 
systems will also exhibit Ehrenfest differences of order the system size on the time-scale 
( 5.16| ). Unfortunately, these authors have not explicitly reported how their results scale 



with increasing system size. 



5.7.1 Liouville Correspondence Does Not Breakdown 

In section I emphasized that the possibility of a breakdown of the correspondence prin- 
ciple was based on demonstrating that macroscopic quantum-classical differences arose on 
an experimentally accessible time-scale. While I have provided evidence above that this is 
indeed the case for the Ehrenfest differences, the magnitude of which grow to the system 
size on the short time-scale ( |5.16| ), I have also demonstrated, in Chapters 3 and 4, that 
the quantum-Liouville differences that arise on this short time-scale do not scale with the 
system size. In particular, in Chapter 3 it was shown that the Liouville break-time grows 
logarithmically with the system size only when applied to differences that remain a factor 
{IT'/·J) smaller than the system size. The presence of such small quantum-classical differ- 
ences in the description of a macroscopic body is not an indication of a breakdown of the 
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correspondence principle, since this principle requires only approximate agreement between 
the quantum and classical theories in the macroscopic limit. 

5.7.2 Can Decoherence Restore Ehrenfest Correspondence? 

In the model that I have considered, each spin is a subsystem of a larger, isolated system. 
The time-evolution of each subsystem may be described, in general, by some non-unitary 
map as a result of the interactions with the other subsystem. The initial coherent state 
for each subsystem rapidly evolves into a nearly random mixture. These properties may be 
considered to arise as a result of the "decoherence" provided from interactions with the other 
subsystem. The important point is that the quantum environment provided by the other 
system only serves to increase the state width of each subsystem, and therefore actually 



leads to a faster breakdown of the Ehrenfest conjecture (5^), but certainly does not help to 
preserve that conjecture. Since I have formulated a statement of the correspondence problem 
specifically in terms of the conjecture ( |5.5| ), as a constraint on the system coordinates, 
these results demonstrate explicitly that decoherence considerations do not circumvent the 
breakdown of Ehrenfest correspondence. 

5.7.3 Consequences for Interpretation of the Quantum State 

As explained in the Introduction, the inapplicability of the conjecture ( p.5| ) to the class of 
chaotic macroscopic bodies, does not entail a failure of the correspondence principle, but a 
failure of a particular interpretation of the correspondence principle. This view is reinforced 
by noting that the conjecture ( |5.5| ) is not implied by Born's postulate, P{q) = \^p{q)\'^, but 
consists of an independent assumption. (However, Born's postulate, combined with the 
assumption (^^), does entail (|5.5|).) 



The Ehrenfest correspondence conditions, (5^) and are motivated by a partic- 

ular interpretation of the quantum formalism, namely, the view that quantum mechanics 
can provide a complete description of an individual physical system. Here it is useful to 
distinguish between two different interpretations of the quantum state [^: 

(i) a pure state provides a complete and exhaustive description of an individual system; and. 



(ii) a pure state provides a description of certain statistical properties of an ensemble of 
similar ly prepared systems. 
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The results of Chapters 3 and 4 provide strong evidence in support of interpretation 
(ii), whereas the argument of the present chapter suggests the possibihty of experimentally 
refuting interpretation (i). 

5.8 Is Quantum Mechanics Complete? 

According to the general argument provided in the Introduction, the quantum expectation 
vàlues may be expected to deviate from the Newtonian coordinates on a time-scale that 
may be short compared to physically relevant time-scales for a wide class of chaotic macro- 
scopic systems. If one is tempted to maintain that the quantum expectation vàlues actually 
describe the dynamics of a single chaotic macroscopic system, as suggested by the Ehrenfest 
Correspondence conjecture (p.5|), then one is lead to the implausible view that significant 
deviations from classical energy conservation should be expected to arise, very rapidly, in 
the case of chaotic macroscopic systems. Since something as fundamental as energy conser- 
vation for an isolated system is likely to be confirmed under a wide variety of experimental 
conditions, it is natural to assume that it is the "expectation value trajectory" predicted 
from the conjecture ( |5.5| ), rather than the Newtonian trajectory, that will be experimentally 
invalidated. 

On the assumption that individual macroscopic systems do not violate energy conserva- 
tion on the time-scale ( |5.16| ), the argument outlined in this chapter leads to the conclusion 
that quantum mechanics is unable to provide a complete description of the motion of in- 
dividual macroscopic systems. Therefore, quantum theory does not provide a complete 
description of reality. This is the conclusion derived by Einstein, Podolsky, and Rosen 
(EPR) a long time ago in a very different way In order to derive their conclusion EPR 



adopted the following principle |34]: 



A necessary condition for a complete theory is that "every element of physical reality 
must have a counterpart in the physical theory. " 

In this chapter I have argued that the coordinates of macroscopic chaotic systems are not 
correctly described by the quantum expectation vàlues since the latter deviate from these 
coordinates in a physically implausible and experimentally significant manner. But since 
this interpretative framework for quantum theory fails to provide a complete description of 
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macroscopic physical reality, is there then any physical interpretation of quantum theory 
that can meet the challenge of providing a complete and correct description of the chaotic 
dynamics of individual macroscopic bodies? 



Appendix A 

Nonclassical Moments for SU(2) 
Coherent States 



Ideally we would like to construct an initial classical density that reproduces all of the 
moments of the initial quantum coherent states. This is possible in a Euclidean phase space, 
in which case all Weyl-ordered moments of the coherent state can be matched exactly by 
the moments of a Gaussian classical distribution. However, below we prové that no classical 
density Pc{0-, (f) that describes an ensemble of spins of fixed length |J| can be constructed with 
marginal distributions that match those of the SU(2) coherent states ( 2.21| ). Specifically, 



we consider the set of distributions on with continuous independent variables G [0, tt] 
and (/) G [0, 27r), measure d/i = smOdOdcj), and subject to the usual normalization, 

dp p,{e,(t>) = i. (A.i) 

2 

For convenience we choose the coherent state to be polarized along the positive z-axis, 
P — b) j)(jij|· This state is axially symmetric: rotations about the z-axis by an arbitrary 
angle (j) leave the state operator invariant. Consequently we require axially symmetry of the 
corresponding classical distribution, 

Pe(e,0)=PcW. (A.2) 

We use the expectation of the quadratic operator, (J^) = j{j + 1), to fix the length of 
the classical spins. 



(J2)c = ViÜ· + l)· (A.3) 
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Furthermore, the coherent state \j,j) is an eigenstate of Jz with moments along the z- 
axis given by (J") = j" for integer n. Therefore we require that the classical distribution 
produces the moments, 

{J2)c=r- (A.4) 

These requirements are satisfied by the 5-function distribution, 

Pv{0) = — ^-5-, (A.5 
zvr sm 60 

where cos^o = Í/|J| defines Oq- This distribution is the famihar vector model of the old 
quantum theory corresponding to the intersection of a cone with the surface of the sphere. 

However, in order to derive an inconsistency between the quantum and classical moments 
we do not need to assume that the classical distribution is given explicitly by ( [A.5D ; we only 
need to make use of the the azimuthal invariance condition ( [A.2| ), the length condition (A.3), 
and the first two even moments of ( |A.4| ). 

First we calculate some of the quantum coherent state moments along the x-axis (or any 
axis orthogonal to z), 

(J™) = forodd m 

(4) = j72 

(4) = 3jV4-i/4. 

In the classical case, these moments are of the form, 

(J™)^ = J dJz J d(l)pc{e)\Jr cos"" {(l))sm"' (6). (A.6) 

For m odd the integral over vanishes, as required for correspondence with the odd quantum 
moments. For m even we can evaluate ( [A.6| ) by expressing the r.h.s. as a linear combination 
of the z-axis moments ( [A.4| ) of equal and lower order. For m = 2 this requires substituting 
sin^(^) = 1 — cos^(0) into ( |A.6D and then integrating over (p to obtain 

(4')c = 7T J dJzp,{e)\3\^ - 7T J d JzPcieM^ cos\e) 

= |J|/2-(j2)/2. 



Since (J|) is determined by (A_^) and the length is fixed from ( [A.3| ) we can deduce the 
classical value without knowing p{0), 

{Jl)c = j/2. (A.7) 
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This agrees with the value of corresponding quantum moment. For m = 4, however, by a 
similar procedure we deduce 

{4)c = 3iV8, (A.8) 
that differs from the quantum moment ( J^) by the factor, 

S4 = 1(4) - {jX\ = |3iV8 - j/M, (A.9) 

concluding our proof that no classical distribution on can reproduce the quantum mo- 
ments. 



Appendix B 

Wigner-Weyl Representation of 
Quantum Mechanics 



The Wigner function is constructed from a quantum state operator p through the transform 




pM,P, = ^ / - y/2|p(í)k + y/2)e*^P/^ (B.l) 

The Wigner function, pw{q,p,t), which is a real "c-number" function, plays the role of 
the quantum state operator in the Wigner-Weyl representation of quantum mechanics. It 
generally takes on negative vàlues in non-vanishing regions of its "phase space" domain, 
{q,p} G R, and therefore may not be interpreted as a probability distribution. 

Each Hermitian operator, A{q,p), in the Standard representation, is also associated with 
a "c-number" function, 

AUq,P) = J dy{q - y/2\A{q,p)\q + y/2)e'^ï'/^ (B.2) 

Observable quantities are obtained from the prescription, 

{A{q,p)) = J dq J dp pu,iq,p) Ay,{q,p). (B.3) 

The Wigner-Weyl representation is an exact representation of quantum mechanics in the 
sense that, 

J dq J dp pyj{q,p) Aw{q,p) = Tr[A{q,p)p], (B.4) 
for any Weyl-ordered operator A{q,p). 
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The time-dependence of the Wigner function can be expressed in the following form |71], 



dpw{q,p) 

dt 



2 



2n 



(-1)" 92n+l^^(^^^)^2n+l^ 



(B.5) 



with an obvious connection to the time-dependence of the Liouville ensemble. Here {•, •} 
denotes the classical Poisson bracket. As noted by Takahashi [^] and others, the quantum 
corrections to the Liouville flow (sometimes called the "Moyal terms" ) do not trivially vanish 
as "/i — > 0" , in spite of their explícit proportionality to powers of h. The Wigner distribution 
has an implícit dependence on h, and factors of (q/h)"' may be expected from the d^/dp^ 
derivatives in Eq. ( B.5 ). 

A clearer picture of the origin of quantum corrections to the Liouville picture is obtained 
by considering the differences that arise in the time-development of "observable" quanti- 
ties. The time-development of the quantum expectation value {A{q,p)), corresponding to a 

after 



time-independent Weyl-ordered operator, A{q,p), follows from Eqs. (B.5) and (| 
integrating by parts, and assuming the surface terms vanish (for physically realistic states), 



d''pw{q,p)d^A 

Qpk Qpj 



0, 



(B.6) 



it follows that, 

d{A^,{q,p)) 
dt 



2 



2n 



(_1^'T·+1 S^n+l^ ^2n+l^ 



{H^,A^} + ^^i^^j (2n + l)! 5p2"+i 5g2"+i " 



(B.7) 



This expression demonstrates that quantum dynamical corrections to the classical Liouville 
flow do not arise unless the Hamiltonian contains at least cúbic powers of g, and, further- 
more, explícit corrections do not arise unless the observable A contains at least cúbic powers 
of p. These conclusions were determined previously by Ballentine and McRae in the specific 



case of the Henon-Heiles model |11]. As before, the quantum corrections may not be ex- 
pected to vanish trivially in the limit "7i — > 0" since the Moyal terms may contain implícit 
dependence on ti. 



Bibliography 



[1] H. Ammann, R. Gray, I. Shvarchuck, and N. Christensen, Phys. Rev. Lett. 80, 4111 
(1998). 

[2] V.I. Arnold, Mathematical Methods of Classical Mechanics (Springer-Verlag, New York, 
1978). 

[3] L.E. Ballentine, Rev. Mod. Phys. 42, 358 (1970). 
[4] L. E. Ballentine, Phys. Rev. A 44, 4126 (1991). 
[5] L. E. Ballentine, Phys. Rev. A 44, 4133 (1991). 
[6] L. E. Ballentine, Phys. Rev. A 47, 2592 (1993). 

[7] L.E. Ballentine, Y. Yang, and J.P. Zibin, Phys. Rev. A 50, 2854 (1994). 

[8] L.E. Ballentine in Fundamental Problems in Quantum Physics, Eds. M. Ferrero and A. 
van der Merwe (Kluwer Acadèmic Publishers, 1995). 

[9] L.E. Ballentine, Quantum Mechanics (World Scientific, Singapore, 1996). 

[10] L.E. Ballentine and J. Zibin, Phys. Rev. A 54, 3813 (1996). 

[11] L.E. Ballentine and S.M. McRae, Phys. Rev. A 58, 1799 (1998). 

[12] L.E. Ballentine, Phys. Rev. A 63, 024101 (2001). 

[13] L.E. Ballentine, private communication. 

[14] J.E. Bayfield and P.M. Koch, Phys. Rev. Lett. 33, 258 (1974). 
[15] G.P. Berman and G.M. Zaslavsky, Physica 91 A, 450 (1978). 

104 



BIBLIOGRAPHY 



105 



[16] M.V. Berry, N.L. Balasz, M. Tabor, and A. Voros, Ann. Phys. 122, 26 (1979). 

[17] M.V. Berry and N.L. Balasz, J. Phys. A: Math. Gen. 12, 625 (1979). 

[18] L. Bonci, P. Grigolini, Phys. Rev. A 54, 112 (1996). 

[19] F. Borgonovi, I. Guarneri, F.M. Izrailev, Phys. Rev. E 57, 5291 (1998). 

[20] P. Boyd, G. Midhn, R. Gilmore, and H. Solari, The Astrophysical Journal 431, 425-431 
(1994). 

[21] A. Braun, P. Gerwinski, F. Haake, H. Schomerus, Z. Phys. B 100, 115 (1996). 

[22] G. Casati and B.V. Chirikov, in Quantum Chaos: Between Order and Disorder (Cam- 
bridge University Press, Cambridge, 1995). 

[23] G. Casati and B.V. Chirikov, Phys. Rev. Letters 75, 350 (1995). 

[24] B.V. Chirikov, F.M. Izrailev, and D.L. Shepelyansky, Sov. Sci. Rev. C2, 209 (1981). 

[25] B.V. Chirikov, F.M. Israilev, and D.L. Shepelyansky, Physica D 33, 77 (1988). 

[26] J.L Cirac and P. ZoUer, Phys Rev Lett 74, 4091 (1995). 

[27] D. Cory and Y. Weinstein, private communication. 

[28] S.R. Degroot, La Transformation de Weyl et la Fonction de Wigner: une Forme Alter- 
native de la Mecanique Quantique (Les Presses Universitaires de Montreal, Montreal, 
1974). 

[29] D.P. DiVincenzo, Phys. Rev. A 51, 1015 (1995). 

[30] J.R. Dorfman, An Introduction to Chaos in NonEquilibrium Statistical Mechanics 
(Cambridge University Press, Cambridge, 1999). 

[31] B. Eckhardt, S. Fishman, J. Keating, O. Agam, J. Main, and K. Muller, Phys. Rev. E 
52, 5893 (1995). 

[32] A.R. Edmonds, Angular Momentum in Quantum Mechanics (Princeton University 
Press, Princeton, 1957). 

[33] P. Ehrenfest, Z. Physik 45, 455 (1927). 



BIBLIOGRAPHY 



106 



[34] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 777 (1935). 

[35] J. Emerson and L.E. Ballentine, Phys. Rev. A 63, 052103 (2001); quant-ph/0011020. 

[36] J. Emerson and L.E. Ballentine, Phys. Rev. E 64, 026217 (2001); quant-ph/0103050. 

[37] A. Farini, S. Boccaletti, and F.T. Arecchi, Phys. Rev. E. 53, 4447 (1996). 

[38] M. Feingold and A. Peres, Physica 9D 433 (1983). 

[39] M. Feingold, N. Moiseyev, and A. Peres, Phys. Rev. A 30, 509 (1984). 

[40] M. Feingold, N. Moiseyev, and A. Peres , Chem. Phys. Lett. 

[41] M. Feingold and A. Peres, Phys. Rev. A 34, 591 (1986). 117, 344 (1985). 

[42] J. Ford, G. Mantica, G. Ristow, Physica D 50, 493-520, (1991). 

[43] J. Ford and M. Ilg, Phys. Rev. A 45, 6165, 1992. 

[44] J. Ford and G. Mantica, Am. J. Phys 60 (12), 1086, (1992). 

[45] R.F. Fox and B.L. Lan, Phys. Rev. A 41, 2952 (1990). 

[46] R.F. Fox and T.C. Elston, Phys. Rev. E 50, 2553 (1994). 

[47] R.F. Fox and T.C. Elston, Phys. Rev. E 49, 3683 (1994). 

[48] H. Frahm and H.J. Mikeska, Z. Phys. B 60, 117 (1985). 

[49] P. Gaspard, Chaos, Scattering and Statistical Mechanics (Cambridge University Press, 
Cambridge, 1998). 

[50] D. Giulini, E. Joos, C. Kiefer, J. Kupsch, I. Stamatescu, and H. D. Zeh, Decoherence 
and the Appearance of a Classical World in Quantum Theory (Springer, Berlin, 1996). 

[51] H. Goldstein, Classical Mechanics, 2nd ed. (Addison-Wesley, Reading, 1980). 

[52] J. Gong and P. Brumer, Phys. Rev. E 60, 1643 (1999). 

[53] M.C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer- Verlag, New 
York, 1990). 



BIBLIOGRAPHY 



107 



[54] F. Haake, M. Kus and R. Scharf, Z. Phys. B 65, 361 (1987). 

[55] F. Haake, Quantum Signatures of Chaos (Springer-Verlag, New York, 1991). 

[56] S. Habib, K. Shizume and W.H. Zurek, Phys. Rev. Lett. 80, 4361 (1998). 

[57] B.S. Helmkamp and D.A. Browne, Phys. Rev. E 49, 1831 (1994). 

[58] B.S. Helmkamp and D.A. Browne, Phys. Rev. Lett. 76, 3691 (1996). 

[59] Ph. Jacquod and J.-P. Amiet, J. Phys. A 30, 2963 (1997). 

[60] R.V. Jensen, Nature 355, 311 (1992). 

[61] E. Jocs and H.D. Zeh, Condensed Matter 59, 223-243 (1985). 

[62] A.R. Kolovsky, Phys. Rev. Letters 76, 340 (1996). 

[63] B.L. Lan and R.F. Fox, Phys. Rev. A 43, 646 (1991). 

[64] J.G. Leopold and I.C. Percival, Phys. Rev. Lett. 41, 944 (1978). 

[65] R.L. Liboff, Introductory Quantum Mechanics (Addison-Wesley, Reading, 1989). 

[66] A.J. Lichtenberg and M.A. Lieberman, Regular and Chaotic Motion (Springer-Verlag, 
New York, 1992). 

[67] A. Messiah, Quantum Mechanics (Wiley, New York, 1966). 

[68] E. Merzbacher, Quantum Mechanics (Wiley, New York, 1970). 

[69] G.J. Milburn, S. Schneider, and D.F.V. James, lon Trap Quantum Computing with 
Warm Ions, Fortschritte der Physik, 801-810, 48 (2000); quant-ph/9908037. 

[70] P. Miller and S. Sarkar, Phys. Rev. E 60, 1542 (1999). 

[71] J.E. Moyal, Proc. Cambridge Phys. Soc. 45, 99 (1949). 

[72] A. Perelomov, Generalized Coherent States and Their Applications, (Springer-Verlag, 
New York , 1986). 

[73] A. Peres, Phys. Rev. A 30, 504 (1984). 



BIBLIOGRAPHY 



108 



[74] A. Peres, Phys. Rev. A 30, 1610 (1984). 

[75] A. Peres, Quantum Theory: Concepts and Methods (Kluwer, Dordrecht, 1993). 

[76] A. Peres, Phys. Rev. E 53, 4524 (1996). 

[77] M. Price, S. Somaroo, A. Dunlop, T. Havel, and D. Cory, Phys. Rev. A 60, 2777 (1999). 

[78] D. T. Robb and L. E. Reichl, Phys. Rev. E 57, 2458 (1998). 

[79] R. RoncagUa, L. Bonci, B.J. West, and P. GrigoUni, Phys. Rev. E 51, 5524 (1995). 

[80] M.E. Rose, Elementary Theory of Angular Momentum (Wiley, New York, 1957). 

[81] J.J. Sakurai, Modern Quantum Mechanics (Benjamin-Cummings, Menlo Park Cahf., 
1985). 

[82] R. Schack and C. Caves, Phys. Rev. Lett. 69, 3413 (1992). 

[83] R. Schack and C. Caves, Phys. Rev. Lett. 71, 525 (1993). 

[84] R. Schack, G. D'Ariano, and C. Caves, Phys. Rev. E 50, 972(1994). 

[85] R. Schack and C. Caves, Phys. Rev. E 53, 3257 (1996). 

[86] R. Schack and C. Caves, quant-ph/9810050. 

[87] R. Schack, Phys. Rev. A 57, (1998); quant-ph/9705016. 

[88] L.I. SchifF, Quantum Mec/tamcs(McGraw-Hill, New York, 1968). 

[89] A.I. Schnirelman, in KAM Theory and Semiclassical Approximations to Eigenfunctions 
(Spinger-Verlag, New York, 1993). 

[90] M. Tabor, Chaos and Integrabüity in Nonlinear Dynamics (Wiley, New York, 1989). 

[91] K. Takahashi, Prog. Theor. Phys. Supp. 98 (1989). 

[92] S. Tomsovic and E.J. Heher, Phys. Rev. E 47, 282 (1993). 

[93] E.P. Wigner, Phys. Rev. 40, 749 (1932). 

[94] J. Wisdom, S. Peale, and F. Mignard, Icarus 58, 137 (1984) 



BIBLIOGRAPHY 

[95] G.M. Zaslavsky, Phy. Rep. 80, 157 (1981). 

[96] W.H. Zurek and J.P. Paz, Phys. Rev. Lett. 72, 2508 (1994). 

[97] W.H. Zurek and J.P. Paz, Phys. Rev. Letters 75, 351 (1995). 

[98] W.H. Zurek and J.P. Paz, Physica D 83, 300 (1995). 

[99] W.H. Zurek, Physica Scripta T76, 186 (1998). 



